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Stress and Deformation in Conical Vibrating Plate’ 


Kazuo IKEGAY At 


Distributions of stress and deformation in a conical thin shell with axial applied force 
are obtained analytically, and the stress and deformation are found to be concentrated on 
the inner and outer edges of the shell. The axial stiffness of a conical plate is extremely 
large compared with that of a flat plate; therefore, forming by press working has proven 
to be a very useful method of producing diaphragms which vibrate as a rigid piston. 
The relations between stress, axial stiffness, and slope of a conical shell are also considered 


in this paper. 


1. Introduction 


As mentioned in a previous paper,‘’? it is 
thought that the axial stiffness of the conical 
part of telephone transmitter and receiver 
diaphragms is large enough for diaphragm 
to vibrate as a rigid piston within the desired 
frequency range. 

But the magnitude of the stiffness depends 
on the slope of the conical shell; and the 
resonant frequency of the bending vibration 
of the conical shell, which may determine 
the upper limit of useful frequency range of 
the diaphragm, depends on the magnitude of 
the stiffness. And especially to prevent creep- 
ing of telephone receiver diaphragms, it is 
desirable that the static stress and strain 
produced by the magnetic force should not 


* MS in Japanese received by the Electrical Communi- 
cation Laboratory on 27 April, 1960. Originally 
published in the Kenkya Zituyéka Hokoku (Electrical 
Communication Laboratory Technical Journal), N.T.T., 
Vol. 9, No.6, June pp. 681-689, 1960. 

+ Station Apparatus Research Section. 


be concentrated on one any part of the coni- 
cal shell; and furthermore the magnitude of 
the static stress and strain should be as small 
as possible. 

Therefore. the internal stress and strain of 
the conical part of diaphragms produced by 
static axial force have been considered, the 
stress on a conical thin shell has been ana- 
lyzed statically, and the distributions of 
bending moments and normal components of 
displacements which may be the most im- 
portant factors for a conical diaphragm were 
obtained. The relations between these magni- 
tudes and the slope of the conical shell were 
considered, and its axial stiffness was also 
investigated. 


2. Fundamental Equations and Their 
Solutions 


Theoretical studies on the thin shell have 
already been carried out by many authors. 
(2)~(6) Tn this paper, the results for the axial 
symmetric case were applied to the conical 


thin shell. 
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As in the papers by the authors referred 
to above, we shall use the following notation 
for the element of the conical shell shown 
in Fig. 1: 


T,=meridian membrane force per unit 
length of shell section 
N,=shearing force in plane perpendicular 
to meridian plane per unit length of 
shell section 
M,=bending moment in meridian plane 
per unit length of parallel circle 
w=normal component of displacement 
1=distance along cone element, measured 
from vertex 
a=angle between axis of rotation and 
an element of cone 
h=thickness of shell 
E=Young’s modulus 
v=Poisson’s ratio. 


(e) 


Fig. 1—Element of conical shell. 


In general the vertical angle 2a of the 
conical diaphragm is nearly 180°, therefore 
in studying fatigue, creeping, and axial stiff- 
ness, it is sufficient to take into consideration 
only the bending moment M, as stress and 
the normal component w as displacement. 
Accordingly the attention is payed to M, and 


w as stress and strain throughout this paper. 
According to the references’”~°? it is 
known that the following equations hold: 


if= —N/Itana+C tan Qa, C1) 


Eh dw 


> 4 Co dw 
12(01—»?) \ dP? 


i rie 
dw __tan’a ee dU _ U 5) (3) 
dl ~~—« Eh dl? dl I : 

where C is a constant which is determined 

by the conditions of the external force, and 


Putting 
K=DA a T (5) 
and ° 
og 1201—»?) 
ype 
h?tan?a ’ (6) 


the function U satisfies the following ordinary 
differential equation of the fourth order. 


fU.. 2: @U 29 (Us S2ae a 
dit | x dB 2 de! Fn 
(7) 
The general solution of Eq. (7) is: 
U=NiI=A,( +-bers—— bei’ x) 
x 
; 2 
+ Af —beix— — ber’ x) 
x 
: Z 
+ Ad( —keix— — ker’ x) 
x 
+ Ad —kerx+ = kei’ x) 
x 
(8) 


where A, As, As, and A, are coefficients 
determined by the boundary conditions, and 


ber x+j bei x=Jo(aj 47) = h(x J 47) 
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ker x+j kei x=Ky(4/ +7) 


ber’ x= 4 ber EG etc. 


Outer Edge 


Fig. 2—Conical shell. 


3. Bending Moment and Displacement 


The coefficients A’s in the general solution 
U may be determined by the boundary con- 
ditions as follows: Assume the conical shell 
is as shown in Fig. 2, the distances from the 
vertex to the inner and outer edges are J, 
and /;, and axial forces p are applied on 
these edges; while the slope at the edges is 
maintained stable by clamping the edges. 
Then the boundary conditions are 


i y) ; 
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N=} at /=/, and h, @) 

dw 

Wi at —Jn and lo, (10) 
and 

T=4— cota at I=), and i (dip 


therefore using Eq. (1), it follows that 


Cae cosec? a@. (12) 
Wigs 


Substituting Eq. (8) into Eq. (3), 


dw #tana [ ners ee as i ae \ 
i? ae), | Ai( bei x - ber *y 


+ Ad( —ber «+ 2 bei’ x) 
% 

+ Ad{ —Kerx+ eae, x) 
x 


+Ad kei x+ & ker’ )+C] 
ii % 
(13) 
Using the results that Eqs. (8) and (13) satis- 
fy the boundary conditions Eqs. (9), (10), and 


(12), the coefficients A’s can be determined, 
hence we have 


Me 
M=A, + bei’ x4 ee —bei x— ce ber’) | 
x x MG 
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+A, eer’ x ss a—»( +berx— bei’«) | 
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xt sin?a@ 2x’ 
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and integrating Eq. (13) and eliminating its 
integration constant, we have 


tan” a / 
GA Ee exibersc— ber») 
w= ORh [ 1 


+ A;(—<x bei’ 1+2bei x) 


+ A3(—x kei’ x+2kei x) 


2p log | 
— / —— ——< . 
+ A,(—x ker’ «+2ker x) + Usaha 


(15) 


Putting w;=w)i=n, and w.=w)i=n, the stiffness 
for the axially applied force shown in Fig. 2 
may be considered nearly equal to p/(w2—w). 
This is easily obtained from Eq. (15). 


4. Numerical Examples 


In this section numerical examples are given 
for the telephone diaphragms made of 0.06 
mm thick duralumin which are shown in 
Lekbeaeee 


Fig. 3—Conical vibrating plate. 


In Fig. 4 and Fig.5 the distributions of the 
bending moment M, and the normal com- 
ponent of displacement w along the cone 
meridians are shown for r=7.5mm., H=3 
mime and) for 7=2.5.mm., 
spectively. 

From these curves it may be seen that 
the bending moment ™, concentrates on both 
inner and outer edges, and at the same time 
its magnitude is very small at the middle 
part. This tendency is evident especially for 
the long meridian of cone J,—1, as shown 


Series 


Inner Edge Outer Edge 


Fig. 4—Distributions of M; and w 
G=7. 5mm —3 mm): 
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Fig. 5—Distributions of M, and w 
Cars sino, J6=S}iva), 


in Fig.5. For comparison the curve of M, 
for the short meridian (r=10 mm., H=3mm.) 
is also given in Fig.6. The reversals of 
bending moments shown in these curves are 
caused by clamping the edges. The bending 
moment at the inner edge is larger than at 
the outer edge; therefore it is desireable to 
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consider the magnitude of the bending moment 
at the inner edge in discussing the magnitude 
of the stress in a cone. These considerations 
on the distribution of the bending moment 
are true also for the normal component of 
displacement, but it may be seen that the 
strain at the inner edge is smaller than at 
the outer edge. 


=O 


21Me / P 


—0,02 
Bil 
—0.04} 


—0.05 


Inner Edge Outer Edge 


Fig. 6—Distribution of M, (v=10 mm, 
H=3mm). 


Fig. 7—Conical vibrating plate with bent 
meridian. 


Next, the effect of the bent meridian shown 
in Fig.7 has been investigated. It has been 
assumed that the radius of curvature of the 
meridian (taken as 50mm _ here) is much 
larger than the other dimensions or that the 
meridian is almost straight. Then assuming 
that the same differential equation, Eq. (7), 
also holds in this case, and putting the 
angles 2a, and 2a, as in Fig. (7) (2a=a,+ as), 
then the boundary conditions for such a 
shell reduce to 


_ p sina, ee 
a ree te | 
(16) 
_ p sina ie 
- Qn sina i=l, 
and 
Ey COS ‘ea 
a 2x cosa lah | 
(17) 
Sy op Fe | 
~ Qn cosa ee 


instead of Eqs. (9) and (11) respectively. 
Therefore using boundary conditions (16) and 
(17), 


a PD COS(G Oy) Ps cosa e (18) 


~ De sin2a Qn sin2a 


and another boundary condition (10), the 
coefficients A’s and therefore the bending 
moment M, can be obtained. 

The dotted line along the full line labelled 
M, in Fig. 4 indicates the distribution of the 
bending moment in this case (r=7.5mm, H 
=3 mm, radius of the curvature of meridian 
=50 mm), which shows a slight decrease at 
the inner edge and for slight increase at the 
outer edge. 

There is another point of view in the in- 
vestigation of these numerical results. In 
Fig.3 holding the height H constant and 
varying the inner radius 7 (or angle of slope 
@) the bending moments at both edges and 
the axial stiffness behave as shown in Fig. 8. 
From these results it is noticed that the 
moment at the outer edge varies slightly but 
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Fig. 8—Bending moments and axial stiffness vs. 


inner radius of conical shell. 
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Fig. 9—Bending moments and axial stiffness vs. 


height of conical shell. 
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that of the inner edge increases greatly with 
the slope becomes gentle, and the rate of 
increase is smaller than that of the inverse 
of inner radius 1/r. 

Then let us vary the height and slope of 
conical shell holding the inner radius constant. 
The bending moments at the edges and the 
axial stiffness for 7r=7.5mm are shown in 
Fig. 9. In this case, the bending moments 
increase slightly but the axial stiffness de- 
creases with the slope becoming to be gentle. 
For zero height shell or flat plate“??? 


p 4x E fh 1 
= 2 , 
We—W, 3 1-0? @ j= ee pb)? 
Me: 


G9) 


and this value indicated in Fig. 9 is extremely 
small in comparison with those of moderate 
height or slope. This shows that the axial 
stiffness of flat plate can considerably be in- 
creased by only press working even if its 
curvature is slight. 


5. Conclusion 


Distribution of bending moment and dis- 
placement in conical thin shell applied the 
external axial force has been investigated ana- 
lytically and numerically, and these magnitudes 
are found to clearly concentrate on the inner 
and outer edges of the shell. The axial stiffness 
is obtained by the difference of displacements 
at both edges, and its magnitude is much 
larger than that of a flat plate. 

Then it has been proved the well-known 
experience that the slight swelling out the 
plate may be effective in keeping its vibration 
pattern like that of a rigid piston in the de- 
sired frequency range. 

These results may also give‘useful data for 
designing conical telephone diaphragms. 
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Some Investigations on the Conical Vibrating 
Plates for Telephone’ 


Kazuo IKEGAYAt and Masayuki MURAKAMIt 


For the purpose of understanding the properties of conical telephone diaphragms, larger 
conical vibrating plates for experiments were made and their characteristic constants were 
measured. From these results the “Frequency constant” and “Stiffness constant have been 
defined as very useful constants for the design of vibrating plates. Their values were de- 
termined experimentally. Then by observations of modes of vibration, it has been concluded 
that the radial corrugations on the conical vibrating plates are unneccessary. 

In the next place, the vibrating plates for new telephone transmitter and receiver have 
been investigated and the effectiveness of the frequency constant and the stiffness constant 
above mentioned have been verified. And it has been known that the method of fixing the 
magnetic armature into the vibrating plate for receiver has influence upon tts modes of 
vibration or area factor, then a most useful method has been newly proposed for fixing 


the armature. 


1. Introduction 


Metallic diaphragms pressed conically are 
used in telephone transmitters and receivers. 
Because of this process of press working it 
is thought that the central part of a dia- 
phragm acts mechanically as a rigid additive 
mass and acoustically as a sound radiator in 
contact with the medium. But in designing 
a conical diaphragm it is necessary to deter- 
mine the fundamental resonant frequency of 
the vibrating system which consists of a rigid 
additive mass and a flat ring plate. 

In this paper, some investigations on the 
conical diaphragms of the No.4 type tele- 
phone are described, and to facilitate more 
reasonable design procedures the “Frequency 
constant” and the “Stiffness constant” which 
will be very useful have been defined. Alumi- 
num conical vibrating plates twice the size 


MS in Japanese received by the Electrical Communi- 
cation Laboratory on 6 April, 1959. Originally published 
in the Kenkyt Zituyéka Hokoku (Electrical Communi- 
cation Laboratory Technical Journal), N.T.T., Vol. 9, 
No.5, May pp. 431-443, 1960. 

Station Apparatus Research Section. 


+ 


of the actual transmitter and receiver dia- 
phragms used in the No.4 type telephone 
were produced for convenience in measuring 
modes of vibration. 

The vibrating-displacement frequency re- 
sponse and the frequency constant were both 
measured in a vacuum using the double-size 
vibrating plates. Also, the stiffness constants 
of the flat ring plate have been determined 
experimentally. These constants have been 
compared with the constants of the dia- 
phragms actually used in telephones. 

By observing the modes of vibration of 
the experimental vibrating plates, we deter- 
mined that radial corrugations on the conical 
section give no beneficial effects; furthermore, 
a new useful method of fixing the magnetic 
armature to the diaphragm of the new tele- 
phone receiver was developed. 


2. Frequency Constant 


As mentioned above, to enable easier and 
more accurate measurements to be made, 
double-size conical vibrating plates were 
fabricated and their properties were investi- 
gated. Fig.1 illustrates these double-size 
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vibrating plates. The plates are made of 0.1 
mm thick aluminum, and thirty different types 
with various values of D, d and H, were pre- 
pared. Moreover for comparison, several ad- 
ditional plates of the same dimensions but 
with seven radial corrugations were fabricat- 
ed. The excitation was electromagnetic to 
a small permendur disk which was pasted 
on the center of each plate as a magnetic 
armature. 


, 10 
« 1004 - a 


« 108] ¢ 


Fig. 1—Experimental vibrating plate. 


The measurements of the vibrating-displace- 
ment frequency response of the double-size 
plates were all made in vacuum because eli- 
mination of the effects of the radiation im- 
pedance to the acoustic medium was desired. 

Fig.2 shows the typical displacement 
frequency response at the center of the vibrat- 
ing plates. Next, the effective mass and the 
fundamental resonant frequency of these 
plates were obtained by the additive-mass 


for No.16 only 


Relative Displacement . 


01 025033 


OSTO7 
Frequency (kc) 


Fig. 2—Displacement frequency responses of 
experimental vibrating plates. 


method. It goes without saying that in this 
case the magnetic gap between the exciting 
pole and the armature was extended enough 
to eliminate the influence of the negative 
stiffness in the gap. As the added mass a 
small brass disk was pasted on the center of the 
plate. From these measurements it was shown 
that the reciprocal square of the fundamental 
resonant frequency is directly proportional to 
the additive mass within the region where 
the weight of the added mass does not ex- 
ceed several times the weight of the vibrating 
plate. In other words, the modes of vibrat- 
ion of the plate are not transformed even if 
an added mass of several times that of the 
plate is added to it. Therefore it follows that 
the fundamental resonant frequency even 
when an added mass is used can be easily 
calculated. Fig. 3 shows an example of these 
data. 


500 


= 488 c/s 


= =e = 4 4 
2 il 0 1 2 3 4 5 6 wf 
Added Mass (g) 


Fig. 3—Determination of effective mass 
by additive mass method. 


From these results we may realize a conical 
vibrating plate with any desired resonant fre- 
quency. In general, for a clamped circular 
plate of radius a and thickness ¢, the natural 
frequencies are 

Ogre tA) 


oS Te eee oo 


fmn : natural frequency of (m,n) mode 
Ann: eigenvalue of (m,n) mode 

E : Young’s modulus 

o : density 

go : Poisson’s ratio 


It would be very convenient if the same 
expression could be applied to the funda- 
mental resonant frequency of the conical 
vibrating plate under consideration. For this 
reason, it 1s assumed the fundamental resonant 
frequency can be expressed by‘! 

Ee lae 
De V pd=at) 2 


which is similar to Eq. (1). Then from ‘this 
equation it is possible to design a conical 
vibrating plate with the desired fundamental 
resonant frequency by a value of the constant 
a determined experimentally. The authors 
call this non-dimensional constant @, the 
“Frequency constant” of a conical vibrating 
plate. The results which were obtained with 
every vibrating plate show that the funda- 
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mental resonant frequency fo of a_ vibrating 
plate is dependent not upon the diameter d 
of the central part of the vibrating plate or 
its height H, but upon the inner diameter D 
of the vibrating ring. The relation between 
the Frequency constant a and the ratio of 
the inner and outer diameters of the ring, 
u=D/Do, is given in Fig. 4, in which the 
dotted lines indicate the 90% _ confidence 
zone. It is shown that the values of this 
constant for the double-size vibrating plates 
agrees with the values for the diaphragm of 
telephone receivers and transmitters. 


+ 


++ = - 
New Transmitter Dilaphragm — a 
4 e 6 


0) 02 04 06 08 10 
H=D/Oo 


Fig. 4—Frequency constant of telephone 
diaphragms. 


3. Stiffness Constant 


In general a preformed vibrating plate is 
considered as a vibrating ring with a rigid 
part at its center (see Fig. 5). When the outer 
and inner diameter of the vibrating ring are 
2a and 2b respectively, and its thickness ¢, 
the axial static stiffness of this system is 
given by ®2@ 


We JU al 


S=— ee 
3 l—-o @& 


Ay? a! 
1-y?— nea Clog )* 
(3) 
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where 


p=b/a 
FE: Young’s modulus 
o : Poisson’s ratio. 


Vibrating Ring 
\ 


\ / 
\ te 
NC Hy 


7/Rigid Part (Sound Radiator ) 


Fig. 5—Model of preformed vibrating plate. 


On the other hand, we may obtain the ef- 
fective stiffness s) of a vibrating plate from 
Fig 3, and, assuming that s) is expressed by 


ee, eB 
= Bo ot Dy e 


similarly to Eq. (3). The authors called this 
non-dimensional constant {) the “Stiffness 
constant” of a conical vibrating plate, whose 
value depends not ugon d and H but on D, 
and therefore on v. The relation between 
f) and # is given in Fig. 6, where the dotted 
lines indicate the 90% confidence zone. The 
values corresponding to f, in Eq. (3) are 


16z 1 
2 
3 1-/2 --*, Clog 12)” 
—# 


aes so (5) 
aH ne — ae a ae aistaterete 
30 30 


x=1—p=(a—b)/a, 


and these values are shown in the same 
figure. The difference between the experi- 
mental values of 8 and the calculated ones 
from Eq. (5) seems to be due to the pure 
statical investigation of Eq.(5) or to round- 
ness of the shoulders indicated by the arrows 
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Fig. 6—Stiffness constant of telephone 
diaphragms. 


in Fig.1. It is also observed that values of 
the constants of the double-size experimental 
diaphragms agree well with the values of 
telephone transmitter and _ receiver dia- 
phragms. 


4. Modes of Vibration 


In order to conjecture the characteristics 
and especially the area factors of a preform- 
ed vibrating plate, observation of the modes 
of vibration of the plate is quite fast; further- 
more this observation is necessary. 

Figs. 7 and 8 show examples of the modes 
of vibration of vibrating plates at their funda- 
mental resonant frequencies. The former is 
for a plate without corrugations and the 
latter for a plate with radial corrugations on 
the conical parts, but both of the figures are 
quite alike and show that the conical part of 
a plate vibrates almost as a rigid piston. 
Figs.9 and 10 also show examples of two 
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plates with and without radial corrugations at 
the second resonant frequency. These modes 
seem to be identical, similar to the case of 


Figs. 7 and 8. 
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Phase Angle 


Fig. 7—Modes of vibration of experimental 
vibrating plate (D=90, d=35, H= 


12.0, without radial corrugations). 
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Fig. 8—Modes of vibration of experimental 
vibrating plates (D=80, d=35, H= 
9.0, with radial corrugations). 
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Fig. 9—Modes of vibration of experimental 
vibrating plate (D=70, d=15, H= 


8.0, without radial corrugations). 
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Fig. 10—Modes of vibration of experimental 
vibrating plate (D=80, d=35, H= 
7.5, with radial corrugations). 
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Table 1 
DIAPHGRAGMS FOR NEW TELEPHONE RECEIVERS AND THEIR CONSTANS 


©)Same dimensions as type © 


t=0.075 Epicote resin applied 


. mean measured value 3386 
’ calculated value 3454 


No.| fo mo | So. A fo | So. mnllenre No.| fo 70 Siena leks 
1 | 1710 | 0.66 |0.761x104 0.64 | 1 | 1564 | 1.03 ]0.994x108}0.72] 4 | 1495] 1.08 |0.952x10% 0.77 
2 | 1707 | 0.69 0693x104 0.75 | 2 | 1554] 1.03 0991x104 0.74] 5 | 1522] 106 |0.968x104] 0.77 
3 | 1726 | 0.67 10.787X10% 0.67 | 3 | 1555} 1.04 j0.992x104 0.73] 6 | 1560| 1.06 |1.07 x 109 0.78 
fo, mean measured value 3498 .Mean measured value 2925 : 
rae calculated value. 3700 | “calculated value’ 3454 | ig meet 2701 | 
(A T7 i 7 7 
©) Lat < Somddmeisions 25 Hoe © ©) 1458 s + €) Same dimensions as type © 
y Ke} 
\ | fo. eet ee Raa 
aa ——— nal A 
picote resin applied += 0.07 Epicote resin applied 
: foe 3% K_|No.| 7 | 7% 50 ZNO EF | a SS 
1872 | 0.72 |0.995x10°} 0.71 1867. 0.82 | 1.13x108| 0.79 | 1 | 1868) 0.81 ]1.11x108) 0.75 | 4 | 1897] 0.82 )1.16x10* | 0.79 
1789 | 0.73 |0.921x10°| 0.72 1819 | 0.80 | 1.04x10°| 0.80 | 2 | 1904] 0.72 |1.03x10®) 0.77 | 5 Seinen ee 
1878 |0.76|1.06 x10} 0.71 3 | 1853] 0.86 |1.16X10? | 0.81 


-mean measured value 3169 
“calculated value 3454 


fo 


+ 


Notes. fo’: fundamental resonant frequency (*%), fo: fundamental resonant frequency of plate only (sv), 


t: thickness (mm), mp: effective mass (g), 


From these results we could not perceive 
any effect of radial corrugations on _ the 
conical part of vibrating plate. So it can be 
concluded that so far as the slope of conical 
part of vibrating plate is within this limit 
bending vibration will be prevented. 


5. Diaphragms for the New Telephone 
Receiver 


On the basis of the above investigations 
five duralumin diaphragms for the new 
telephone receiver, which are illustrated in 
Table 1, have been made. In this table the 
measured values of fundamental frequency, 
effective mass, effective stiffness, area factors, 
and the value of f) calculated from the fre- 
quency constant a are also given. 


5.1. Vibrating Displacement and Sound 
Pressure 


Examples of vibrating displacement at the 
center of these diaphragms in vacuum and 
sound pressure in a 7.7 c.c. coupler are shown 


So: effective stiffness (dyn/cm), K: area factor. 


in Figs.11 and 12. These two curves are 
similar to each other, so we may expect that 
the diaphragms under test are useful for the 
new telephone receiver whose desired fre- 
quency range is upper limit about 4000 cycles 
per second. 


5.2. Modes of Vibration at the Funda- 
mental Resonant Frequency 


As the modes of vibration of a diaphragm 
have an effect upon its effective area; name- 
ly, on the sensitivity of the receiver, we have 
observed the modes of vibration at fundamental 
resonant frequency in order to obtain the data 
to judge its quality. In Figs.13 to 17 ex- 
amples of modes of vibration are shown, and 
it can be seen that these modes are de- 
pendent upon the existance of radial corru- 
gations on the conical part of the diaphragm. 
For diaphragms without corrugations (Figs. 
13(a), 16 and 17) the displacement at the 
circumference of the flat central part or at 
the part fixing the armature is remarkable, 
and for plates with corrugations (Figs. 13(b), 
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In Vacuum 


Relative Displacement 


if; 2 S 5 7 
Frequency (\kc 


Fig. 11—Vibrating displacement at the center 
of new telephone receiver diaphragm 


(Type ©). 


14 and 15) the magnetic armature is fixed 
extremely rigidly by the corrugating process 
increasing the stiffness sufficiently to prevent 
deformation. From these investigations we 
see that if we can fix the armature without 
decrease of stiffness, the modes of vibration 
will be much improved, and therefore increase 
in the area factor and improvement of the 
sensitivity of the telephone receiver may be 
expected. But it must be noticed that the 
radial corrugations were found to be useless 
in the former investigations on the experi- 
mental vibrating plates. This is because in 
the former case the magnetic armatures were 
pasted on these plates, and so no decrease 


F THE ELECTRICAL COMMUNICATION LABORATORY 


| 7 Fexe-Coupler 


Relative Sound Pressure 


Frequency (kc) 


Fig.12—Sound pressure produced by new 
telephone receiver diaphragm (Type 


©). 
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(a) Type @ (b) Type ® 


Without Radial Corrugation With Radial Corrugations 


Fig. 13—Modes: of vibration of new telephone 
receiver diaphragms. 


of stiffness was caused by the process of fix- 
ing the armature. Accordingly the modes of 
vibration were ideal. On the other hand, the 
process of fixing the armature to the dia- 
phragm of a telephone receiver causes a de- 
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Fig. 14 Fig. 15 


Fig. 14—Modes of vibration of new telephone 
receiver diaphragm (Type ©, with 
ring armature and radial corrugations). 


Fig. 15—Modes of vibration .of R-4 receiver 
diaphragm (with radial corrugations). 
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Fig. 15 Fig. 16 


Fig. 16—Modes of vibration of new telephone 
receiver diaphragm (Type ©). 


Fig. 17—Modes of vibration of new telephone 
receiver diaphragm (Type ®). 


crease of stiffness, but this decrease may be 
prevented by the use of corrugations. That 
is to say the contribution of corrugating the 
conical part of the diaphragm is not to stiffen 
it but to prevent a decrease of stiffness. 


5.3. Method of Attaching the Armature 


As mentioned above, the deformation of 
the circumference of the flat central part of 
telephone receiver diaphragm without radial 
corrugations is extremely large. Therefore in 
order to recheck this deformation, it was 
necessary to stiffen the deformed part by 
some suitable way and then to observe the 
modes of vibration. In the experiments this 
part was stiffened by applying Epicote resin. 
The fundamental resonant frequency, effective 
mass, effective stiffness and area factor 
of vibrating plates with resin are given in 
Table 1. Example of the modes of vibrat- 
ions measured for these plates with resin, at 
the fundamental resonant frequency, are also 
shown in Figs.18 and 19. As the constants 
and the modes of vibration can be varied 
and the area factor raised by applying resin, 
a suitable method of fixing the armature by 
which no decrease of stiffness is caused must 
be obtained. Fig. 20 shows a new method 
of fixing the armature which resembles the 
method of attaching caps to beer bottles. 

The experimental data of the constants of 
telephone receiver diaphragms for which the 
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Fig. 18 Fig. 19 


Fig. 18—Modes of vibration of new telephone 
receiver diaphragm (Type ®’, Epicote 
resin applied, see Fig. 16). 

Fig. 19—Modes of vibration of new telephone 
receiver diaphragm (Type ©’, Epicote 
resin applied, see Fig. 17). 
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Table 2 


CONSTANTS OF DIAPHRAGMS FOR TELE- 
PHONE RECEIVER WITH NEW METHOD 
OF FIXING THE ARMATURE 


No. t So’ Mo So K 


1 ONOGSELOZONmE Ont 0.802108 80.9 
ORCL LO ZI ORG i ON/0Z< 102 Deo 


Oo bw 


0.065 1745 0.79 0.949x10% 80.9 
4 0.063 1615 0.79 0.813108 80.2 
5 0.061 1546 0.75 0; WOW Ge 7). I 


fo’: fundamental resonant frequency (¢s), 
t: thickness (mm), 

mo: effective mass (g), 

Sy : effective stiffness (dyn/cm), 

K : area factor. 
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Fig. 20—Diaphragm for new telephone re- 
ceiver with new method of fixiug 
the armature. 


new method of a fixing the armature was 
adopted are given in Table 2, and it may be 
seen that this new method is as effective in 
increasing the area factor as the application 
of resin is. This method is adopted in the 
new telephone receiver. 


5.4. Creeping 


For the electromagnetic receiver the creep- 
ing of the material of the diaphragm caused 
by magnetic attraction is a serious problem 
because it is required that there be no static 
variation of the air gap length. From this 
reason the magnitude of creeping of the dia- 
phragm (see Fig. 20) for the new telephone 
receiver to which the new method of fixing 
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the armature applied was measured. The 
applied load was 1 kg weight at three points, 
but the magnitude of creeping was saturated 
at under one micron so it can be concluded 
that the creeping in this case is negligible 
(see Fig. 21). 
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Fig. 21—Magnitude of creeping of the dia- 
phragm for the new _ telephone 
receiver with new method of fixing 
the armature. 


Table 3 


DIAPHRAGMS FOR NEW TELEPHONE TRANS- 
MITTERS AND THEIR CONSTANTS 


"TSS Se 
32.2 ¢ | : 
35 - | 48 53 a 
a Ag Ve) pe Le wo 
t=0.052 |With Corrugations t=0.051 
fo =745 fo = 697 
My = 0.195 my = 0.196 
Sy CS NO? SS OOn 
A =719 =e 
© 469-4 © es 46 6——— 
meee -o=32 2 — | 
| ese | eee 1 | & A8.58- a 
With Corrugations +=0051 t= 0.062 
—| 
I SHS} fo = 890 
M = 0.188 My = 0.226 
Sep CIAO S<il(oy So = 7.0810 
kK =/05 jG fils) 
Notes. 
fo: fundamental resonant frequency of diaphragm 
only (3), 


t: thickness (mm), 
my: effective mass (g), 

Sy: effective stiffness (dyne/cm), 
K: area factor. 
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6. Diaphragm for Telephone Trans- 
mitter 


Four duralumin diaphragms for the new 
telephone transmitter, which are shown. in 
Table 3, have been made. In this table the 
mean values of the fundamental resonant fre- 
quency, effective mass, effective stiffness, and 
area factor are given for groups of three to 
five plates. 


In Vacuum 


Relative Displacement 


OLS am! 2S 5 10 
Frequency (kc) 


Ol 02703 


Fig. 22—Vibrating displacement at the centers 
of diaphragms for new _ telephone 
transmitters. 


6.1. Vibrating Displacement 


An example of the vibrating displacement 
at the center of these diaphragms by electro- 
magnetic excitation in vacuum are shown in 
Fig. 22, and from this result we may expect 
that the diaphragms under test are useful for 


the new telephone transmitter desired fre- 
quency range upper limit is about 4000 cycles 
per second. 
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Fig. 283—Modes of vibrations of new telephone 
transmitters diaphragms. 


6.2. Modes of Vibration at the Funda- 
mental Resonant Frequency 


Some examples of modes of vibration are 
shown in Fig. 23. In this case the ratio of 
the inner and outer diameters y is very small 
as compared with the ratio w of the receiver 
diaphragms (see Table 3), and consequently 
the conical part of the diaphragm acts almost 
as a rigid piston and its bending vibration is 
produced only at the flat ring part. The 
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modes of vibration of a diaphragm with an 
electrode of 0.08 g weight is shown in Fig. 
24, and those for a diaphragm with the 
electrode and the carbon granule in Fig. 25. 
They were excited by the standard receiver 
RM-303 through a small coupler. 
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Fig. 24 Fig. 25 


Fig. 24—Modes of vibration of new telephone 
transmitter diaphragm with movable 
electrode (Type ©). 


Fig. 25—Modes of vibration of new telephone 
transmitter diaphragm with movable 
electrode and carbon granules (Type 


©). 


7. Conclusion 


In order to better understand the properties 
of conical diaphragms for telephone trans- 
mitters and receivers, some experimental 
double-size conical vibrating plates have been 
made, and their fundamental resonant fre- 
quency, effective mass, effective stiffness, etc. 
have been measured. From these results the 
“Frequency constant” and the “Stiffness con- 
stant”? which are very useful for the design 
of telephone diaphragms have been defined, 
and their values were experimentally deter- 
mined. Then by observing the modes of 
vibration it has been concluded that radial 
corrugations on the conical diaphragm are 
unnecessary. 

Next investigations on ' diaphragms for 
telephone transmitters and receivers have been 
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conducted and asa result the effectiveness 
of the frequency constant and the stiffness 
constant have been verified. And it is shown 
that the method of fixing the magnetic 
armature to the telephone receiver diaphragms 
has an influence upon the modes of vibration 
and area factors, and an extremely good 
method of a fixing the armature has been 
proposed and its effectiveness has also been 
confirmed. 

Of course there are many other important 
problems concerning telephone transmitter 
and receiver diaphragms other than the pro- 
blems discussed above; namely the second 
resonant frequency and the modes of vibrat- 
ion at its frequency, etc. Moreover, there are 
many unascertained aspects about the press 
working; the minimum amount of working 
and their order to obtain the required 
accuracy, the magnitude of stroke of one 
press working, the design of press die, etc., 
and these questions call for further investi- 
gations. 
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3. Study of Auditory Confusion 


When the response mode of the listeners 
based upon the “identification” method of 
signal-quality judgment is adopted as our 
principal method of quality measurement, the 
most intensive interest of our studies is 
focused upon the relation of correspondence 
between an ensemble of timbre-signals as 
stimuli and an ensemble of timbre-qualities 
as responses, that is, on the problem of the 
“auditory confusion.” For example, if we 
want to understand more concretely the 
quality characteristics of BED so far ex- 
plained, we must come to the points to be 
decided: Under a particular band-cutting 
distortion, what is the degree of confusion of 
a particular voice signal with which this 
signai corresponds to quality responses? Also, 
how does this type of confusion correspond 
to the pattern deformation due to this dis- 
tortion? This does not mean a_ one-sided 
study of the listener’s response based only 
upon “correct” judgement as found in the 
past tranditional articulation studies, but takes 
into consideration also the minute character- 
istics appearing in the “incorrect”? judgment. 
We’can thus attain novel and clear interpre- 
tations of the “quality” itself and even the 
“distortion.” The term “quality modification” 
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+ Nagoya University. 
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which has been used corresponding to the 
pattern deformation stems from the special 
study of the auditory confusion presented 
here. 


3.1. Confusion in the Distortionless 
Condition 


When clarifying the quality phenomena 
from the viewpoint of correspondence between 
a deformed signal pattern and the listener’s 
quality response, there are some problems 
about the fact that signal and quality do not 
show one-to-one correspondence even under 
the distortionless condition, where no physical 
distortion exists. Although the present study 
is aimed at a quality study of speech under 
physical distortion, this subjective quality 
study in the distortionless condition will be 
considered as the basis of our study, because 
it is thought to indicate directly the concrete 
mechanism of quality formation. The quality 
to be dealt with first is voice quality. 


3.1.1. Confusion Matrix and Voice 
Quality Configuration 


The confusion matrix is obtained by tabu- 
lating the frequency of correspondence ob- 
served in each combination of signal and 
quality. Table 2 shows the confusion matrix 
for BED of experiment No.6. The capital 
letters vertically aligned on the left side of 
the confusion matrix are signal symbols 
showing kinds of voice at the sending point. 
The lower-case letters horizontally aligned in 
the under part are quality symbols showing 
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Table 2 


VOCAL QUARITY CONFUSION MATRIX IN THE 

ALL-Pass CONDITION (EXPERIMENT NO. 6) 

a Ne Se ee eee eee 
All-pass Condition 


| 
~ | t eh || aa h Di 
S| =e =a! 
| 
al 362 116 16 6 | 500 
Si 84 322 92 iz, | 500 
M’ 21 72 Sy) Oe | 500 
H 4 19 3 464 500 
es 
SS 471 529 476 524 2000 
Notes: 


Capital letters=Transmitting signals. 
Lower-case letters=Receiving qualities. 
Blocks with heavy outline show definite sub- 
groups. 
Blocks with dotted outline show somewhat 
ambiguous subgroups. 
el=elision. 
Voices=All male voices. 


listeners’ response at the receiving point. 
The figure in cells shows the aggregate of 
frequency of correspondence between signals 
and qualities for four listeners. It is seen 
from the confusion matrix that there are 
some trends in confusion occurrence and fur- 
ther that any two voices which are adjacent to 
each other in such arrangement are mutually 
related most closely. This arrangement of 
voices is called “voice-quality configuration.” 
It can be shown that the voice ensemble 
adopted here is further divided into 3. sub- 
groups as indicated by the smaller blocks in 
the table. When we examine the arbitrary 
combination of any two voices, it is possible 
to consider the “reciprocity” of confusion. 
For example, as seen from the figures in the 
table, the trend of confusion from “T” to 
“s’”, is more marked than that from “S’” to 
“t?, and in such a case the confusion is said 
to be nonreciprocal. For a combination of 
“S’” and “M’”’, on the other hand, there is 
fairly similar frequency of correspondence, 
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between “S’’?>“m’ and’ “‘M'’—s" 5 andmit 
may be said that the reciprocity holds to a 
certain extent. On the contrary, lack of re- 
ciprocity between “M’”’ and “H” is most 
marked. From the fact that the reciprocity 
does not hold in some combinations, it follows 
that there seems to be a certain kind of bias 
in judgment of listeners regarding criteria of 
the differentiation of two voices of similar 
voice quality. How such a bias is establish- 
ed can not be immediately answered, but it 
may be inferred that the reason for the lack 
of reciprocity that seems to take place be- 
tween “M’” and “H’”, for example, is that a 
clue to quality judgment of “H” is better 
grasped than that of “M’”’. In the case of 
“T? and “S'¥ “Tas only similar’ tor ac 
quality of “s’”, while “S’’ has qualities similar 
to both “t” and“m’”. “It “can then Be ex 
plained that the confusion frequency of “S’” 
is shared by two different confusion-objects, 
thus resulting in nonreciprocal confusion. 
In any case, the manner of configuration by ° 
which voice qualities are arranged in the 
ensemble reveals, so to speak, some matter 
and substance regarding the listeners’ per- 
ception of voices and consequently this is 
intimately related to the listeners’ subjectivity. 
To begin with the discussion of reciprocity 
and lack of reciprocity in confusion so far 
described is just to stress how the idea of 
“direction” is indispensably important in the 
analysis of confusion phenomenon. For further 
details of fundamental interpretation of con- 
fusion, reference must be made to the Lite- 
rature (38). 


3.1.2. Relation between the Envelope 
Shape of the Vocal Pattern and 
Voice-Quality Arrangement 


Now we pass into the comparative obser- 
vation between voice-confusion in the dis- 
tortionless condition and the vocal pattern of 
individual voices. [Experiment No. 6.) 

Here analytical study was done by utilizing 
the sonagraph. Fig.9 is obtained by averag- 
ing five vowel patterns in each voice given 
by the sonagraph, in which, however, charac- 
teristics of higher frequency emphasis of the 
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Fig. 9—Base envelopes found in each vocal- 
pattern (Experiment No. 6). 


analizer are duly compensated for and adjust- 
ed to exactly flat characteristics. 

The almost linearly decreasing curves are 
drawn as shown in the pattern figures, by 
tracing the basis of the principal glens of 
the pattern as one of the simple ways of 
characterizing the shape of these patterns. 
These curves thus obtained can not be said 
to correspond directly and immediately to 
some physical concrete substance. However, 
they may be considered as bearing upon 
some direct expression of vocal pattern charac- 
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teristics which is governed at least by the 
vibration modes of the vocal cord and _pro- 
bably further by the resonance position and 
selectivity mode of an articulatory organ. 
Fig. 10(A) shows an averaged pattern of 
each vocal pattern, and Fig. 10 (B) gives only 
the envelopes obtained after discarding all 
peak structures from each vocal pattern. We 


propose that these envelopes be designated 
the “base envelope.” 
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Fig. 10—Comparison of averaged form of vocal 
pattern and “base envelope” with regard 
to four individual voices. 


When comparing each voice, it is well to 
notice that the slopes of the envelope curves 
run-in’ the order of 1’; S.-M and hie 
and further that this order conforms to the 
voice-quality arrangement in the confusion 
matrix as already explained. From this fact 
it is assumed that the voice-quality classificat- 
ion has something to do with such a simple 
extraction of shape characteristics of vocal 
patterns. Careful observation naturally can 
be carried out by comparatively examining 
the shape of the envelopes of the four indi- 
vidual vocal patterns, instead of their base 
envelopes, and the shape of an averaged 
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pattern of four voices. For example, the 
shapes of the patterns of “S’” and “M”’ are 
very similar to the averaged one, while the 
other two, “T” and “H”, particularly “H”, 
are not similar to the averaged one. Joint 
consideration of this fact and the ambiguity 
of voice judgment observed in the confusion 
matrix yields the conclusion that the voices 
having the patterns similar to averaged one 
are most ambiguous in regard to voice-quality 
judgment. From this it follows immediately 
that when listeners try to classify voices in 
a given ensemble for memorizing them in 
their mind, the principal scheme of their 
voice-classification does not depend on a 
particularly localized frequency band of pat- 
terns, but relies possibly upon silhouette 
shapes outlined throughout the whole patterns. 
The extractive representation of pattern by 
base envelope as shown in Fig. 10(B) is 
simple and effective, but can not be consider- 
ed to furnish us sufficiently full information, 
although it is even of interest to note that 
rough arrangement of voice quality is possible 
very reasonably in this way. 
Lexperiment No. 1) 

The voice ensemble in this experiment con- 


Table 3 


VOICE QUARITY CONFUSION MATRIX IN THE 
ALL-PASS CONDITION (EXPERIMENT NO. 1) 


| f° | 5° | y° | ro | h° ys 
a see ine i a —- ry, 
F° 25209) <3 Cae as 3 || 400 
| Panes Gens = 
S 16 | 321 28 30 5 | 400 
ye 1 24 | 283 34 | 58 | 400 
Cc 1 11 47 | 320 | a 400 
ish 0 | 8 32 | 10 | 350 | 400 
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sists of “F°”’ (male), “S°”’ (male), “C°” (child), 
“v°” (female) and “H°” (female) voices. The 
confusion matrix in the all-pass reference 
condition is shown in Table 3. 

The voice-quality arrangement is shown in 
the order’of “F®1S°, Y°, C’ and’ H° 3) iie 
vocal patterns of the individual voices and 
their averaged pattern are shown in the order 
of the voice quality arrangement in Fig. 11. 
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Fig. 11—Explanation of pattern arrangement 
corresponding to voice configuration 
(Experiment No. 1). : 


The experimental processes to obtain these 
patterns are carried out utilizing a frequency 
analyzer as explained in Chapter 1. The four 
full lines with arrows drawn through these 
pattern figures are: 


VOLUME 8, NUMBERS 11-12, NOVEMBER-DECEMBER, 1960 513 


(a) the line which is formed by the first 
formants of the five vowels and the second 
formants of “A” and “O” (sometimes with 
“Oran 

(b) the line which is formed by the second 
formants of vowels, “I” and “E”. 

(c) the line which is formed by the third 
formants of the five vowels. 

(d) the line which is formed by the fourth 
formants of the five vowels. 

The frequency position of these lines in the 
patterns becomes higher approximately accord- 
ing to the order of voice-quality arrangement, 
with the exception of voice “H°”’, where its 
frequency position of lower formant cannot 
be indicated so definitely as in the cases of 
other voices because of the complexity of the 
pattern shape. It has been found that the 
formant positions approximately coincide with 
the physical stature of the calling subjects 
and seem determinable by the size and volume 
of their articulatory organs. However, the 
position of the higher formants for example, 
the third (c) and the fourth (d) formant lines 
of “H°”’’s voice are relatively high notwith- 
standing her comparatively large physique, 
and these data do not seem connected direct- 
ly with the size of the articulatory organs. 
In voices “C°” and “H®’, pattern parts indicat- 
ed by the (b) line are quite dispersed, re- 
vealing that the positions of the second formant 
of “IT” and “E” are considerably apart from 
each other. This seems to be caused by the 
fact that the calling subjects have pronounc- 
ed “E” with their mouths a little wider open 
than normal. 

Now the relation between vocal-patterns 
and voice-quality arrangement will be ex- 
amined for individual voices. Let us look at 
the outlined forms of pattern envelopes in 
Fig.11. Several facts of importance can be 
seen: First, the form of the intermediate 
voice “Y°” in arrangement is very much like 
that of the averaged pattern; second; “C” 
and “S” are alike; third, the forms of the 
two voices “F°’ and “H°” at both ends of 
the arrangement are the least alike. Similarity 
to the average pattern in these cases is 
noticed not only in the envelope forms as a 
whole, but also in the position of each formant 


as shown by lines (a), (b), (c), (d). These 
present relatively clear and intimate corres- 
pondence between the characteristics of the 
pattern form and modes of grouping in voice- 
quality configuration. As the vocal pattern, 
presently under observation, is too narrow by 
about 25 dB in component amplitude by the 
limitation of oscillographic recording, the de- 
tails of the base envelopes can not be suf- 
ficiently obtained as before. Nevertheless, it 
is of interest, as showing a clue to the me- 
chanims of vocal quality formation, that fairly 
clear correspondence has been found between 
positional relation of the pattern formants 
and voice-quality arrangement. Among various 
voice groups, the “Y°’—““H®” group does not 
apparently correspond to mutual disimilarity 
of these two patterns, but this seems to be 
due to complexity in the lower part of the 
“H®’ pattern, because the “H°’ pattern form 
in its lower part has a part of characteristic 
of “C°” as the smallest physique on one hand, 
and «a. part. of characteristic of Fas) the 
largest physique on the other hand. The 
strong fundamental component of “H°” results 
from the fact that the frequency is already 
situated in the part of the first formants of 
“7? and “U”. Then it follows that in this 
case when pronouncing these vowels, the size 
of the throat cavity is as large as that of a 
male, while the reinforcement of the third 
harmonics may be attributed to the fact that 
when pronouncing other vowels the size of 
the cavity is that of a female or child. The 
voice whose vowel pattern changes from the 
region of a male voice to nearly that of a 
child according to the kind of vowel is pe- 
culiar. The relation between vocal patterns 
and confusion, in such a case, can not easily 
be explained. The separate article (39) may 
be referred to for the details of the charac- 
teristics of voice “H”’. 

[Experiment No. 5] 

The confusion matrix is shown in Table 4. 
The voice qualities are arranged in the order 
Oped. SNe C2. andes. 

Subgrouping by sex, such as “K”’-““M” or 
“C/”_S”, is most conspicuously noticed, then 
follows grouping of “M’-“C”’, which is 
rather nonreciprocal. The respective vocal 
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Tabie 4 


VOICE QUALITY CONFUSION MATRIX IN THE 
ALL-PASS CONDITION EXPERIMENT NO. 5) 
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K’=Male, M=Male, C’=Female, S=Female. 


patterns of individual voices and their averag- 
ed pattern are shown in the order of the 
voice-quality arrangement in Fig. 12. The fre- 
quency analyzer was the same one explained 
in Section 1. 

From the observation of the outlined forms 
of pattern envelopes in this figure, it is easi- 
ly seen, as in the preceding cases, that the 
patterns of voices “C’” and “M’’ having inter- 
mediate voice quality are most similar to an 
averaged pattern. The position and strength 
of formants in the upper part may be con- 
sidered to correspond to the subdivision of 
the ensemble by sex into two groups, “K’’’- 
“M” and “C’”-“S”, as explained in Section 
3. Although the grouping of “C’’-“M” can 
be explained by similarity of pattern forms, 
it is impossible to find any interpretation of 
the nonreciprocal trend of confusion of ‘“C’” 
—>“m” without further minute inspection of 
pattern shapes coming down to the details of 
individual vowels. 

Runing through the data of the above 
three experiments, correspondence between 
the confusion under the distortionless con- 
dition and the relation of vocal pattern 
arrangement is sufficiently traced. It is re- 
cognized from voice-quality arrangement and 
vocal pattern representation that the listeners’ 
classification of voice quality primarily depends 
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Fig. 12—Explanation of pattern arrangement 
corresponding to voice-quality ar- 
rangement (Experiment No.5). 


upon the outlined shapes of pattern envelopes. 
To explain fully, in the cases where the en- 
semble consists only of male voices as in ex- 
periments No.2 and 6, the listeners can 
classify voices only on the basis of the 
simpler shape of the base envelope in vocal 
patterns obtained by excluding peak structures 
from our observation. In experiments No. 1 
and 5 where consciously the sex element 
was brought into the voice ensemble, the 
shapes of the base envelopes do not  corres- 
pond simply to voice-quality arrangements. 
In such cases, however, the listener’s ability 
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to discerne sex difference seems to depend 
primarily upon the position and strength of 
formants in the upper part, and then seem 
to depend upon the relative position of the 
peak structures localized in the frequency 
dimension, as a clue for more delicate quality 
difference. 


3.2. General Treatment of Confusion 
Phenomena 


The effort made in the preceding section, 
is to make clear the matter of voice-quality 
configuration and to trace the mechanism of 
voice-group classification, based exclussively 
upon the confusion matrices in various ex- 
periments. It is by this confusion study of 
voice ensemble that we aimed at deriving an 
interprtation of vocal pattern. However, it 
becomes exceedingly and evidently trouble- 
some to examine the confusion in the most 
detailed phase, because the confusion matrix 
has the same number of degrees of freedom 
as the number of total cells found in its 
matrix. In the preceding section, as an aspect 
of confusion phenomena, the direction of 
confusion is pointed out, and the importance 
of direction for the analysis of confusion 
phenomena is stressed, citing the proposal of 
the concept of reciprocity of confusion. In 
dealing with confusion, other important con- 
cepts are the “amount of frequency” and the 
“equivocation” regarding correspondence be- 
tween signals and qualities. In order to par- 
ticularize the confusion matrix systematically 
and generally, further discussion may be 
made on these fundamental concepts as 
follows. 


3.2.1. Principle of Confusion Analysis 


1) Direction of confusion 

The aim of the analysis of the confusion 
matrix is to clarify the mechanism of quality 
formation by tracing the relation between 
quality response and pattern forms of signals. 
Therefore our main interest regarding con- 
fusion matrix analysis in general distortion 
cases will naturally be focused on the follow- 
ing two questions: (1) To what matter of 


quality responses does the pattern of a_par- 
ticular signal deformed by a certain type of 
distortion correspond? In order words, how 
is a particular signal deformed by distortion 
connected with a set of responses of all 
qualities? (2) For a particular quality, to 
what content of quality response does a_ set 
of patterns of all signals deformed by certain 
distortion correspond. In other words, how 
is a set of signals deformed by distortion 
connected with a particular quality? 
The former question (1) deals principally 
with quality modification of a signal by dis- 
tortion and this problem is necessarily con- 
nected with the concept of “quality deterio- 
ration” or “quality destruction” by distortion. 
The latter question (2) treats exclusively 
quality modification of signals by distortion, 
thus turning out to a particular quality, and 
this problem is necessarily connected with 
the concept of “quality construction” or “quali- 
ty formation“ by distortion. Confusion con- 
sidered from the former standpoint is called 
“outgoing confusion“ which goes from a par- 
ticular signal to an entire group of qualities. 
Confusion considered from the latter stand- 
point is called “incoming confusion“ which 
comes to a particular quality from a group 
of signals. The concept of direction is thus 
introduced by distinguishing between outgoing 
and incoming confusions. These new concepts 
naturally stem from the essential consideration 
of the quality-theory based upon confusion 
analysis. Confusion should not be treated as 
a phenomenon where one signal is simply 
misheard as another. By considering signal 
patterns at the sending point and qualities at 
the receiving point, it becomes possible to 
study the mutual correspondences between 
them along two directions. The quality-theory 
orientated along this- line can have its basis, 
in the interpretation of confusion phenomena. 
2) Frequency of confusion occurrence and 
uncertainty of correspondence 

When analizing the confusion matrix sepa- 
rately in two directions, not only uncertainty 
of correspondence between signals and quali- 
ties, but also frequency of confusion oc- 
currence should be taken into consideration 
as far as we are dealing with an ensemble. 
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When, for instance, a certain distortion in- 
creases not only the uncertainty of corre- 
spondence; but also the frequency of occurrence 
of outgoing confusion from a_ particular 
signal, this distortion should be considered as 
“quality-destructively“ acting on this particular 
signal. But on the other hand, when a cer- 
tain distortion increases not only uncertainty 
of correspondence, but also the incoming con- 
fusion to a particular quality, this distortion 
can be said to be “quality-constructively“ 
acting regarding this particular quality. Intro- 
ducing the concept of confusion direction and 
then evaluating in each direction the occur- 
ring frequency of confusion as well as un- 
certainty of correspondence, it becomes 
possible for the first time to attain both 
generality and simplicity in confusion analysis. 


3.2.2. Interpretation of Phonemic Con- 
fusion Phenomena from the Con- 
cept of Outgoing and Incoming 
Confusions 


It is very important in quality theory to 
introduce the concept of the incoming con- 
fusion as well as outgoing confusion into the 
analysis of confusion phenomena. Confusion 
analysis in the past seems to have been limit- 
ed mostly to what we call outgoing confusion. 
By the combined consideration of incoming 
and outgoing confusions, the problem of. re- 
ciprocity can properly be evaluated. In addit- 
ion the confusion phenomena can be treated 
in more generality than ever, and from these 
basic treatments, several new facts become 
available for the interpretation of quality 
phenomena. For example, the phonemic 
confusion in BED, the quality of which is 
characterized by concentration of quality 
density. The outgoing-confusion characteristics 
can fully account for the process of quality 
destruction, as the nature of the confusion 
on this direction may be considered to be 
strongly related to the eliminated part of the 
pattern of phoneme. On the other hand, the 
incoming confusion may be considered as 
directly bearing upon the residual quality, as 
it has much to do with the remaining part 
of the pattern. It is in this respect that the 
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incoming confusion phenomena are full of 
interest. From these points of view, Fig. 13 
is made to express the importance domain of 
each phonemic quality of five vowels. In the 
figure, phonemic quality characteristics in the 
ordinary sense are fully utilized, being com- 
bined with the data of phonemic confusion 
characteristics: namely, for each vowel, the 
percentage value of correct articulation and 
the ratio of incoming confusion to particular 
phonemes are used jointly in the same 
figure. 
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Fig. 13—Representation of quality distribution 
of 5 Japanese vowels by confusion 
characteristics and quality loss charac- 
teristics (Experiment 2). 
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Even from the results of this experiment, 
where the cutoff conditions were quite rough 
and could not cover the entire frequency 
range required, it is easily seen that an in- 
coming confusion to a particular quality, is 
fairly stabilized regardless of the kind of 
source from which it originates, and further 
that the mode of incoming confusion appears 
very concentrated on the frequency scale. 
It will also be distinctly seen in this figure 
how the position of the quality-balancing 
point relates to the real position of quality 
accumulation. 

The new concept of “negative importance” 
is duly discovered by taking the incoming 
confusion into consideration in addition to 
the outgoing confusion. 
‘Interpretation of the frequency band can _ be, 
from the viewpoint of quality distribution, 
expected to be more delicate and fine. The 
traditional thinking of importance as_ being 
only “positive”, seems to stem from the point 
of view which is limited to outgoing con- 
fusion. Where a given freauency-band position 
corresponds to a quality accumulating point 
of a particular phoneme quality, the confusion 
incoming to this phoneme suddenly increases, 
because the position may be regarded as 
“positively important for this phoneme. On 
the contrary, when some particular band re- 
gion other than the one considered above 
corresponds to a very conspicuous quality 
accumulating point of some phoneme other 
than the phoneme considered, the outgoing 
confusion from this phoneme suddenly in- 
creases on that band region; this means that 
the position of this kind will prove itself not 
only unnecessary but even harmful to the 
phoneme in consideration. The concrete ex- 
amples of the description will be found in 
Fig. 13. “Negative importance,” or the con- 
cept of-“‘nuisance,’ comes out in this way. 
The problem treated above is limited to the 
consideration of phonemic quality only. 
However, the really significant problem lies 
in the interaction between voice quality and 
phonemic quality. The transmission theory 
primarily aiming at one-sided transmission of 
intelligibility only and keeping the voice 
quality elements out of sight, necessarily leads 


Then the quality: 


us to an easy thought to regard voice ele- 
ments as only “redundant” for phonemic 
elements. According to us, to do away with 
nuisance does not mean to save the fre- 
quency bands from an economical viewpoint, 
but to express positively the intelligibility. 
This consideration goes out of the sphere of 
faithful transmission based on the traditional 
idea of “fidelity“ in communication. 


3.2.3. Definition of Confusibility 


In order to analyze confusion phenomena 
in two directions, and to represent occurring 
frequency and uncertainty of correspondence 
quantitatively, it is convenient to use the 
same measure of information amounts as 
used in information theory??“?%”, The 
figures in the confusion matrix show occur- 
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MATRIX CONFUSION 
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ring frequency, from which the statistical 
probability (ratio) can immediately be calcu- 
lated. pare 

If in the generalized confusion matrix 1n 
Table 5 (A), statistical probability concerning 
the correspondence between a_ particular 
signal X; and a particular quality x; is shown 
by p(X, x), we obtain: 


Qj 


DX, Xj) a N 


Occurrence probabilities of this signal and 
this quality are respectively 


5 N.; 
pXD= a y Pay) aN 


For conditional probability between signal X; 
and quality x;, we have 
Qij 


paX=_, PXiMD=N 


from which is obtained the following relation: 
PX, 4) =P KD P@IXD) =P) PK). 


Therefore, as shown before in Table 5(B), 
frequency ratio in the matrix can be com- 
pletely replaced by the statistical probability. 
Using the values thus transformed, the quanti- 
ty to express the uncertainty of correspond- 
ence in the direction from a single to a 


quality group, that is the outgoing confusion 
will be: 


Clea aa LPH ii) log p(x Xi) 


which, further weighted by the probability of 
occurrence, leads to the expression: 


Cou =— LP(Ks, Xi) log p(xjXi) 


The individual outgoing confusibility is de- 
fined by this quantity in which are taken 
into account at the same time the uncertain- 
ty of correspondence and frequency of corre- 
spondence. In the same way, regarding the 
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incoming confusion. 


(Ou Ses 2 DX, x; logp( Xi x;) 


is obtained for a quality x,, and it defines 
the individual incoming confusibility. These 
quantities are measured in the unit of “bit” 
if the base of the logarithms is taken as 2. 
Of course, the sums of individual confusibilit- 
ies of all signals and of all qualities are equal 
to “uncertainty” H(xX) and “equivocation” 
H(X'\x) respectively. The total sums, there 
denoted by Cour and Cin, are called ensemble 
outgoing confusibility and ensemble incoming 
confusibility and calculated by the following 
formulas: 


Cour= DCout, Cin= Chin. 
i j 
Further, the information amount transmitted, 


R, is calculated from the confusion matrix in 
the following manner: i.e., 


R=H(X)—-HACX'x) =H) -—H'X) 
= HX) —Cin= H (x) — Cour. 
To make the occurrence frequency of each 
signal equal and to avoid elision in perception, 
as has been done in our experiments, the 
following relation 


H(X)2 HAH) 


always exists. 
lation 


Thus, we have finally a re- 


Cire Cou: 


The individual confusibility defined above 
generally takes into account both the fre- 
quency and the dispersion of correspondence 
between signal and quality. Strictly speaking, 
however, further consideration concerning the 
occurrence mode of each signal is needed to 
combine the dispersion and the frequency of 
correspondence into a factor directly reflecting 
the quality effect of distortion. But, as the 
occurrence frequency of signals is kept con- 
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stant throughout all of the experiments, the 
outgoing confusibility in our case is reduced 
to the representation of only the dispersion 
of correspondence. Therefore the increase or 
decrease in outgoing confusibility can be con- 
siderd as reflecting simply the degree of 
quality destruction due to distortion. Although 
the incoming confusibility takes into account 
both the dispersion and frequency of corres- 
pondence, it can be also conceived that the 
confusibility in this direction truly represents 
the degree of quality construction. The large 
value of the incoming confusibility means 
that to a considerable extent all signals tend 
to have a particular quality. For the exact 
evaluation of quality effect of each condition 
of distortion, however, it is necessary to 
compare the value of the outgoing confusibili- 
ty with that of the incoming one. When, for 
instance, quality judgment becomes sufficient- 
ly ambiguous as a whole by certain distortion, 
it is quite conceivable that the outgoing con- 
fusibility as well as the imcoming confusibili- 
ty increases. If, for certain distortion, the re- 
lation: 


Cant > Cra 


is obtained, this distortion is considered to 
have a destructive action of a quality x, with 
regard to the signal group; and if the relation 


Chim oO ont 


is obtained, the distortion is considered to 
have a constructive action of the same quality 
regarding the same signal group. 


3.2.4. Representation of Experimental 
Results by Confusibility 


Making use of the data of measurement in 
BED of experiment No.2 as an example, 
the confusibilities defined above are calculat- 
ed from the confusion matrix for each con- 
dition of distortion. The explanation of the 
experimental results based upon this manipu- 
lation of data will now be described. The 
confusion matrices are shown in Tables 6 
and 7 for phonemic quality and voice quality 


respectively. 

(1) Phonemic confusion 

Individual confusibilities for five vowel 
phonemes are calculated in each cutoff con- 
dition and shown in Fig. 14. The results for 
LC and HC distortions are presented in Figs. 
14(B) and (C) respectively. In Fig. 14(B) is 
shown five vowel patterns, each of which is 
averaged for four voices. It will be noticed 
from the general form of confusion curves 
that confusibilities in both directions are very 
different from each other in the particular 
cutting points, and also the  confusibility 
curves take very characteristic shapes accord- 
ing to the kinds of vowel sound. In the 
succeeding paragraphs, we shall explain the 
shapes of curves for each combination of 
vowel sounds and quality. 

Oo fe a 

The characteristic curve in HC distortion 
is more remarked than that in LC distortion. 
Outgoing confusion has rather simple charac- 
teristics in both distortion directions. This 
means that diminuition in phonemic quality 
is not remarkable in spite of the deformation 
of patterns within the distortion range.  In- 
coming confusion in LC distortion becomes 
gradually conspicuous from LC 1.0ke and 
reaches about 1.7 bits at LC 2.5kc. As the 
conceivable maximum value of individual 
confusibility is 2.3 bits, it can be said that 
about 83% of vowels in signal ensemble 
have the quality “i” at this cutting point. 
The formant components between 2-3 kc 
found in each pattern of the Fig. (B) have 
undergone a sort of emphasis influence, as 
the result of sudden cutting of the components 
in the lower band and level recovering of the 
remaining part at the same time. Thus, the 
band-cutting distortions from 1.5 to 2.5 kc are 
said to bring quality-constructive action of “1” 
with regard to vowel ensemble. With regard 
to HC distortion, incoming confusion does not 
surpass outgoing confusion within the range 
of distorting conditions used. As it seems 
too hasty to conclude from this figure that 
neither sudden cutting off in the higher region 
nor emphasis in the lower region contributes 
to forming “i” quality, it is necessary to pre- 
pare more finely and sharply cutting-off steps 
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Table 6 
PHENOME-QUALITY CONFUSION MATRIX IN BED (EXPERIMENT NO. 2) 
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Table 7 


VOICE-QUALITY CONFUSION MATRIX IN BED (EXPERIMENT NO. 2) 


all-pass : Z P Le sit 5 ke 

| y k n f y y k n i wy 
eee aie : 26 ifs aa 162 65 Ay 2 300, 
18 Zl i) i? 300 ee | a 36 i 300 
12 18 250 20 ee 40 | : 5 210 34 te 
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Voice: All male voices. 
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Fig. 14—Representation of phonemic confusion 
in BED and display of vowel patterns 
(Experiment No. 2). 


covering the frequency region lower than 
0.5 ke. The phonemic quality of “I” is duly 
represented by the pattern differences between 
the vowel pattern (in full line) and mean 
vocal pattern (in dotted line) as shown in 
the figure (B), and the differences are observ- 
ed a little concentrated in the frequency re- 
gions of the neighborhood of 2ke and 0.2- 
0.3ke. The components near 2kc, which 
may well be considered as the main formant 
of phoneme “I” are sharply emphasized as 
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Fig. 14—Representation of phonemic confusion 
in BED and display of vowel patterns 
(Experiment No. 2). 


the last residuals of band components in LC 
distortion, while in HC distortion, many lower 
components remain uncut and critical emphasis 
of particular components can not be expected. 

eae 

The confusion characteristics of this vowel 
is similar to those of “I. The strong trend 
in quality construction is recognized from LC 
1.5 kc, and this means, as in the case of “i”, 
that the quality transition to “e’”, i.e., quality 
formation, takes place by emphasizing the 
pattern components of each vowel near this 
cutting-off point (Confusion toward “i” from 
“E” increases at LC 2.5 kc, but incoming to- 
wards “e” decreases there, and thus the 
quality effect of distortion becomes ambiguous). 


In HC distortion, the outgoing confusion is 
generally stronger than the incoming one, 
and particularly, distortion of HC 1.0kce 
seems only to have a destructive action of 
quality “e’. ‘The components near lke as 
shown by the pattern of “E” in the Fig. (B) 
seems to have what we call “negative im- 
portance” of phoneme “E”. To put it con- 
cretely, the second formant necessary in 
quality formation of “e” is cut off and the 
timbre of residual patterns comes to resemble 
“Op gtlensaesPelC: 

The note regarding “i” 
“e” for HC distortion. 

“A? —“a” 

The confusion characteristic of this vowel 
is that the trend of quality formation becomes 
maximum at 1kc in distortions of both di- 
rections. This means that the partial emphasis 
of components in this region of vowel signals 
is most powerfully contributing to quality 
formation of “a”. Concretely speaking, in- 
coming confusion from “E”, “O”, “U”, etc. 
towards “‘a” is increased in this frequency 
region. 

Marked trend in outgoing confusion at LC 
1.5 ke is mainly due to confusions going to- 
ward e”, “i” and “u’, but at this cutting: 
off point the distortion is absolutely destruc- 
tive for phoneme “A”, being reflected by 
outgoing confusion which tends toward “e”, 
“7, “uu”, etc. Thus, the components in this 
region of vowel “A” are considered to be 
contributing to quality “a” only in the “nega- 
tive’ sense. Outgoing confusion is strength- 
ened at 0.5 kc in HC distortion. As seen in 
the pattern of “A” in the figure (B), this 
corresponds to the fact that the negative im- 
portance of phoneme “A” exists in two parts 
just above and below the frequency band 
tke 

(ee One 

Generally speaking, LC distortion is quality 
destructive and HC distortion is constructive 
for this phonemic quality. Neither the deteri- 
oration in phonemic quality of this vowel, 
nor the quality transition from other vowels 
towards “‘o” will take place until LC 0.6 kc. 
Confusion outgoing from this vowel is a little 
conspicuous at the cutting points of LC 1.0 ke 


will also apply to 
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and up. The pattern components in these 
regions 1.0-1.5 kc, having the action of form- 
ing other phonemic qualities, clearly have 
negative importance for the phonemic quality 
“5”. LC 1.0kc and under are the regions 
where imcoming confusion is dominant —a 
fact due to the conspicuous confusion incom- 
ing from “I” and “E” towards “o”. Probably 
the center of phonemic quality of “o” may 
exist at some points lower than 1 kc. 

SSW 

The characteristics of this vowel are a 
little more complicated than those of others. 
The characteristics seem to be particularly 
well reflected in LC distortion. It has, in 
both directions of LC and HC distortions, 
cutting points yielding the quality formation, 
yet those cutting frequencies do not necess- 
arily conform in both directions. Speaking 
in more detail, emphasis of components by 
LC 1.5ke and the elimination of compo- 
nents by HC 2.1 kc and under are all quality 
constructive for phoneme quality “u’. The 
former corresponds to emphasis of components 
near 1.5kc of “A” and “O” patterns, and 
the latter to the elimination of the second 
formant of “T’ and “E’, respectively. After 
all, both the second formant and the first 
one found in 0.3-0.4kce are contributing to 
quality formation of “u”. This mechanism 
of quality formation is more complicated than 
that observed in other vowels. 

In the above, the trend in phonemic con- 
fusion in the BED with level recovering pro- 
cess has been considered by the aid of in- 
dividual confusibility for each vowel. The 
important fact as found therein is as follows. 
Vowel patterns are, after all, of gradually 
diminuishing type of amplitude with increase 
of frequency domain, component emphasis in 
the neighborhood of the cutting point of a 
filter is more effective in lower-cut than 
higher-cut, and in many cases it is better to 
consider that the main function of higher 
cut-off consists in the elimination of compo- 
nents in the attenuation region of a filter and 
not in the emphasis of the remaining band. 
Therefore, a sharp cutoff filter acts on quality 
“destructively” in some regions and “con- 
structively”, in other regions, and so it is 
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possible to distinguish two kinds of frequency 
bands, namely, the part of “positive im- 
portance” and the part of “negative impor- 
tance’. 

Fig. 15 was obtained by plotting the values 
of ensemble confusibility at each cutting point 
as the total sum of individual confusibilities 
of all phonemes. 


ict 
|Phoneme judgment | 
J eee 


T 


fusibilities (bit) 


Con 


Transmission (%) 


Percentages of Information 


Ensemble 


0.1 O12 OS) (OS OLp a 2 3 5 


Cutott Frequency (kc) 


Fig. 15—Ensemble phonemic confusibility charac- 
teristics in BED (Experiment No. 2). 


The reason C,, is generally not smaller than 
Cour is because of the experimental scheme 
as already stated. Except the LC 2.5kc, HC 
distortion shows the more conspicuous con- 
fusion trend than LC. The equality of out- 
going and incoming confusibilities means that 
no directivity is found in the confusion trend 
in an average sense.* (*If a pair of signals 
is considered, confusion between signals may 
have reciprocity. But when a group consist- 
ing more than two signals are considered, 
regressive confusion trend, such as “I”—“E” 
—"“A”—*T”, may be conceivable. In these 
examples, however, there is found no such 
particular confusion trend). In such a case 
confusion occuring may be considered as be- 
ing at random. The difference between both 
confusibilities is the result of the facts that 
the outgoing confusion from a_ particular 
phoneme decreases, or that the incoming 
confusion to a particular quality increases. 
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In these cases, where the difference between 
confusibilities in observed, the confusion may 
be regarded to have a certain trend instead 
of being at random. In short, confusion trend 
reflects in inequality (Cin>Cou.) in these ex- 
periments. The figure showing the trend just 
mentioned can be a clue to interpret the 
frequency construction of phonemic quality. 
Conspicuous difference between confusibilities 
at LC 2.5kc results from the fact that con- 
fusion towards phoneme quality “i” at this 
cutting point becomes definitive throughout 
the vocal ensemble. LC distortion at this 
cutting point can be interpreted to be re- 
markably quality-constructive. 

With the exception of this cutting point, 
greater difference between both confusibilities 
is more generally found in HC distortion than 
in LC. This means that HC _ distortion 
happens more frequently to behave as quality- 
constructive, and furthermore, that there is 
more negative importance in the lower part 
of the vowel pattern than in the upper part. 
This seems to be caused by the greater con- 
tribution of the second formant to phonemic 
quality, especially in the case of vowels “TI” 
and “E”. As has already been stated, the 
difference between the information amount of 
input signals and ensemble incoming con- 
fusibility, R=H(X)—Cin, shows amount of 
information transmitted. In the figure these 
values at each cutting point are shown in 
percentage for the maximum values in the 
vertical axis on the right side. It is of interest 
to note that about 50% of information 
amounts are transmitted at the crossing point 
of Cin curves. 

2) Vocal confusion 

Contents of confusion matrix obtained from 
voice judgment will be treated by the intro- 
duction of individual confusibility as in the 
case of phonemic confusion. In this case, 
however, some confusions are noticed even 
in the distortionless condition, and further- 
more the confusion occurs in a different way 
according to each voice. In Fig. 16 individual 
confusibility in the reference condition is 
shown for each combination of voice signal 
and voice quality. 

The outcome of such characteristics can not 
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Fig. 16—Comparison of individual outgoing 
and incoming confusibilities fer each 
combination of voice signal and 
voice quality in the distortionless 
condition (Experiment No. 2). 


be clarified immediately, but one fact is that 
the contents of voice quality impressed upon 
listeners depends to a certain extent upon 
the kinds of vowels, as will be later explain- 
ed. It may be said that in the case of the 
voices which are classified by the listener in 
a stable manner regardless of vowel kinds, 
the outgoing confusion is less than the in- 
coming one, while on the contrary, in the 
cases of voices the qualities of which are 
impressed on listener only by a particular 
vowel, incoming confusion is less than the 
outgoing confusion. As an illustration, the 
values of difference Cin—Cou, of individual con- 
fusibilities in two directions are calculated 
and shown in Fig. 17 for each vowel sound. 
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Fig. 17—Comparison of dependency of outgoing 
and incoming confusibility upon kinds 
of vowel in voice judgment in the dis- 
tortionless condition. 
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These confusibilities are obtained from the Vocal Pattern 

confusion matrix which are tabulated for [ ae aa es he ce Uhie oe ee 
each vowel. As will readily be understood As. HC HC HC Hoe 8S 


from the figure, voices “K” and “F” which [ 
generally show more confusion in the incom- 

ing direction give little deviation in con- 

fusibility difference by vowel kinds, on the L 
other hand, voices “N” and “Y” which show 
more confusion in the outgoing direction give 
considerable deviation. From the detailed 
inspection of the confusion matrix, it can be 
said that “N”, the quality of which appears 
most clearly only in the vowels “I” and “U”, 
is very liable to be confused with other voices 
when other vowels than “I” and “U” come 
up for judgment. And “Y”, with the voice 
quality especially appearing in “A” and “O”, 
has greater outgoing confusion in other vowels AY, 
than these two. Anyhow, it should be kept SA 
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Fig. 18—Individual representation of confusibility 
of vocal confusion in BED and display 
of vocal patterns of four voices (Experi- 
ment No. 2). ' 


in mind that the individual confusion due to 
BED start from such a stage in the reference 
condition. In Figs. 18 (A) and (C) are shown 
individual confusibility in each cutoff con- 


Individual Confusibilities (bit) 


wa “N”—« dition of LC and HC distortions with regard 
to each combination of voice signal and voice 
0.2 ; quality. Figs. 18 (B) shows the vocal pattern 
Al xe L Mee of four voices, and to each pattern an averag- 
ed mean pattern is superimposed in dotted 

OS ean line. 
- Poe ap oe The general trend of vocal confusion curves 
el ee is such that the shape of the curve is simpler 
Sle ae than that of phonemic confusion and no cut- 
0.0L. — Naina. 3 ting step is recognized which shows so re- 
cee 0.5 0.7 TT markable a trend toward quality destruction 
Cutoff Frequency (kc) or quality construction as in phonemic quality. 

a o——_—__e Ce a x On 


However, it is true that by using the curves 
in two directions, more abundant information 
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Fig. 18—Individual representation of confusibility 
of vocal confusion in BED and display 
of vocal patterns of four voices (Experi- 
ment No. 2). 


regarding the mechanism of quality establish- 
ment is obtainable than that in the ordinary 
judgment curves already stated. 
Yea ay” 

- For this voice, generally speaking, HC 
distortion is almost always quality-destructive, 
and LC distortion is nearly always quality- 
constructive in frequency band between 0.6—- 
1.0kc. With this voice, together with “K,” 
we can observe the abundance of harmonic 
components as characteristics of their pattern. 
It is therefore conceivable that the elimination 


of these components due to HC distortion 
has caused similarity to the patterns ‘“N” 
and “F” which are originally of attenuating 
type of higher components; this seems to 
correspond to quality destruction of this voice 
in HC destortion. In LC distortion, partial 
enhancement of the pattern parts in the band 
between 0.6-1.0 ke corresponds to the trend 
of quality construction of this voice in this 
band region. Then the positive shape of “Y” 
pattern (in a sence of relativeness to a mean 
pattern) in this band region is taken up as 
one of the quality factors of this voice, but 
this localized factor of voice quality cannot 
be detected by HC distortion, because of its 
general nature of quality destruction. It means 
probably that this quality factor of voice is 
not. able to construct voice quality inde- 
pendently and solely, but dependently call 
upon other factors in other band regions. 
Meanwhile its outgoing confusion remarkably 
increases in 2-3 kc in HC distortion, thereby 
revealing the existance of voice quality ele- 
ment in this band, but this localized factor 
can never be detected in LC distortion. 

a lel ce 

LC distortion is nearly always  quality- 
destructive, and HC distortion constructive in 
one part and destructive in another. Increase 
in outgoing confusion at LC 0.4 kc is mainly 
caused by the increase of the confusion to- 
wards “y.” It is conceivable that the domi- 
nant part of pattern in the narrow region 
up to about 0.4kc of vocal pattern of “K” 
is a constructive factor of voice quality “k.” 
After all, with reference to the preceding 
section, it can be said that the voices “K”’ 
and “Y” are similar with respect to the 
abundance of harmonics in their pattern 
spectrum, and this causes the relatively fixed 
stable grouping of “K”-“Y.” For the dis- 
crimination of voice qualities of the voices in 
the same group, the small difference in vocal 
pattern, in other words, the quality factor 
localized in a particular band as just noticed 
above, seems to be a clue. More confusion 
maintains throughout the LC distortion in 
the outgoing direction than in the incoming, 
and in HC distortion, also the same tendency 
is kept at HC 1kc and lower cutting points. 
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These confusion characteristics can be ex- 
plained by considering certain facts. First, 
notice the dispersed mode of energy distri- 
bution in this voice spectrum. This distri- 
bution mode appears in the peculiar shape of 
the vocal pattern for it has the strong funda- 
mental component, and yet at the same time 
the higher components extend up to the 
highest region of any of the four voices. 
This is probably not the pattern characteristic 
due to vocal force in emission intensity. On 
the other hand, observing the confusion 
relatively increases at HC 4.2 kc. This leads 
to the statement that the quality factors of 
this voice exist near this cutting point or in 
the higher’ band as well as in the extreme 
lower band as mentioned first in this section. 
Thus two distortions HC and LC, in spite of 
their reversed directions, both bring much 
outgoing confusion. 

The increase incoming confusion in HC 
2.1kc is mainly caused by confusion from 
voice “Y.” This is seen in the preceding 
section as increase of outgoing confusion of 
"Yat this cutting’ point. Thus, some of 
voice-quality factors of “K’” may be expected 
to also exist near this cutting point as well. 

IN fe ae 

From the general trend of the curves it 
is seen that LC and HC distortions are 
equally quality-constructive for this voice 
quality. It should be remembered that this 
voice has vocal pattern most similar to the 
averaged one and the outgoing confusibility 
in reference condition has the greatest value 
among the four voices. Thus it can be ex- 
pected that the mode of discrimination for 
this voice quality may not be the same kind 
of differential type as observed in the voice 
group “K—Y,” but may be a different and 
rough type. At LC 0.4ke the incoming 
confusion exceeds the outgoing. This can be 
explained as follows. If one ‘eliminates, on 
the figures, the parts under this cutting point 
from the patterns of each voice, he will 
immediately notice that “Y” and “F” lose 
their characteristics of pattern shape. “K” 
neither loses its characteristics of pattern, nor 
is confused with “n” at this cutting point. 
With regard to “N” itself, it is noticed that 
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it loses similarity to averaged pattern. These 
deformations of patterns can be considered 
to correspond to the decrease of outgoing 
confusion and increase ot incoming confusion 
of this voice at this cutting point. Concerning 
the HC distortion, the increase in incoming 
confusion at HC 4.2kc is mainly caused by 
increase of incoming from “K” and “Y.” 
This seems to relate to the fact that the 
vocal patterns of those voices, by this dis- 
tortion, loss their quality factors localized in 
higher frequency region, and begin to re- 
semble the mean pattern. The larger incoming 
confusion in HC 1ke and 0.5ke is mainly 
caused by increase of confusion which comes 
from “K,” and this simply corresponds to 
the similarity in both shapes in the lower 
part of pattern. 

In short, this voice whose vocal pattern is 
most similar to the mean pattern has no 
voice-quality factor localized in a particular 
frequency band. The increase in incoming 
confusion due to BED comes from the elimi- 
nation of localized quality factors of other 
voices. For this voice itself, ambiguity dis- 
appears from the pattern shape in BED, and 
in some cases, the outgoing confusion de- 
creases. 

Sc eS h22 

Generally speaking, LC distortion is quality- 
destructive and HC distortion, on the contrary, 
is quality-constructive. The confusion, the 
trend of which follows the distortion direction, 
may be said to be systematic, and the cor- 
respondence between the trend and pattern 
characteristics can be comparatively easily 
obtained. The characteristics of vocal pattern 
of this voice is the strongest fundamental 
component and the most remarkable dimi- 
nuition of harmonics among the four voices. 
This is the general characteristic of vocal 
pattern brought about by vocal-intensity of 
pianissimo, and it is readily conceivable that 
low cut or higher emphasis of such vocal 
pattern is to destroy its voice quality. On 
the other hand, higher cut or lower emphasis 
is not a cause to make it lose these pattern 
shape characteristics, but is one to give the 
pattern characteristics of this voice to the 
patterns of other voices. This seems a 
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reasonable explanation of confusion trend in 
this case, that is, LC distortion is quality- 
destructive and HC distortion quality-con- 
structive. 

From the analysis of voice quality confusion 
for four voices as above, it is safe to make 
the following conclusions. The discrimination 
or classification of voice quality generally 
dapends upon outlined shapes of envelope of 
vocal pattern. Concerning each individual, 
however, it can be said that among voices 
forming a particular subgroup by voice 
quality depends on the pattern part localized 
in a particular band. When the voice having 
such localized quality factors subjects to BED, 
the outgoing confusion generally exceeds the 
incoming confusion as the result of disap- 
pearance of quality factors. On the other 
hand, the voice having no localized quality 
factors does not show critical increase in 
outgoing confusion at any cutting step, but 
on the contrary, conspicuous increase of 
incoming confusion is observed as the result 
of the quality destruction of other voices 
having the localized quality factors. Attention 
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Fig. 1g¢—Characteristics of ensemble’ confusibility 
5 of voice quality in BED (Experiment 
No. 2). 


should be paid, however, to the fact that 
localized voice-quality factors recognized by 
analysis of such confusion trend are not 
simultaneously detected in the two distor- 
tions of LC and HC, and further that their 
mutual independency is not so difinite as in 
phonemic quality factors. 

Finally, the ensemble confusion character- 
istics calculated as the sum of all individual 
confusibilities are shown in Fig. 19. 

Compared with the case of phonemic 
quality, it is noticed that the slopes of the 
curves are not so irregular. The same trend 
was also observed before in the judgment 
curves. The difference between outgoing and 
incoming curves is much less than in phone- 
mic quality. Some differences are observed 
only in IC .0.6 ke-and,. HCG 21 ker this 
means that the confusion occurring in voice 
ensemble is random or not systematic, and 
it can be said that BED as a whole is 
neither quality-constructive nor quality-des- 
tructive for voice quality. The difference be- 
tween curves at LC 0.6kc means that the 
incoming confusion towards the voice quality 
“vy” becomes systematic to some extent at 
this cutting point. Another difference between 
curves at HC 2.1 ke shows that this cutting 
point is likely to cause comparatively con- 
spicuous confusion trend of “Y”’—“k.” These 
characteristics of curves both reflect the 
localized quality factors of voice “Y,” and 
mean that, with this voice, the voice quality 
is concentrated on the frequency band to 
somewhat higher degree than with the other 
voices. In this figure, the amount of infor- 
mation transmitted is shown in percentage 
in the vertical line on the right side. The 
crossing point of curves presents 40%; about 
0.6 times the value in the reference condition. 


(to be continued to “Review of E.C.L. 
Vol.9, Nos1-2 issue’) 
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Theory of Solution of Chain Molecules’ 


Hirosi OKAMOTO} 


By using the first order Bethe approximation, a theory of the solution of chain molecules 
is developed. Though it is a well-known fact that the first order Bethe approximation 
and the method of the quasi-chemical equilibrium coincide with each other in the case of 
a lattice gas or of the solution of short chain molecules, the author found that this coin- 
cidence does not always hold for the solution of long chain molecules. His theory can 
better explain the experimental observations on the mixing energy and the phase separation 
than the previous theories do. 


1. Introduction be shown that the results can explain the 
experimental data which the previous theories 
Many years have pzssed since the pioneer- could not. Though it is a well-known fact 
ing work of Flory’s‘” quasi-crystalline model that the first order Bethe approximation and 
of chain polymer solutions was developed. the quasi-chemical equilibrium coincide with 
Miller,“ Orr,“ Guggenheim,‘* and other each other in the case of a lattice gas or in 
investigators also studied a similar problem the case of a solution of short-chain mole- 
and succeeded in explaining the general cules, it was found that this coincidence is 
features of chain polymer solutions. not always true in the case of a solution of 
Dispite the apparent success of these long-chain molecules. 
theories, they are not able to predict accurate- 
ly, the detailed behavior observed experiment- 2. The First Order Bethe Approximation 


ally. For example, the experimental values 
of the concentration at the critical mixing 
point are almost twice as high as those the 
theories predict and the concentration de- 
pendency of the mixing energy is more com- 
plicated than the theories predict.‘” 


A lattice space whose coordination number 
is z is considered. The system consists of 
nm, Q@-type chain molecules and n, §-type 
solvent molecules. The a-type chain molecule 
consists of m elements which occupy 2 con- 


To improve these theories, the concentration tiguous latfice Sites, “while: the: 2-tyge | soln 


dependency of the interaction parameter was molecule is a point molecule occupying ONS 
proposed by Tompa‘” and Huggins.“ Hug- site. The first neighboring pairs of two 


gins also introduced additional concentration cena ee) (8-8), and (a~f) contribute 

dependent terms in the entropy part. But the potential energy of sa, ¢ and s,5 respec- 

the author considers that these methods are UNE 

Pee iicill and aindesieakle: For the case of the first order Bethe appro- 
famines paper the “detailed = conneamnene ximation, all the possible configurations that 

of the polymer will be considered and it will can be: illustrated son) asset @f ep ieetice 


consisting of a central sete and its first 

aS Bee eet cant) coordination shell, provided that the first shell 
Ss receive y the ectrica ommunication La- sites a 1 

boratory on 14, September, 1960. Originally published i. ae oe pe nelgibery 

in the Journal of the Physical Society of Japan, ie considered in detail. The group orthese 

Vol. 15, No. 4, April, pp. 650-657, 1960, z+l1 sites is hereafter called a “local system.” 

High Polymers Research Section. The outer part of the local system affects 


= 
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the detailed configurations of the system in 
the follow.ng way. (Fig. 1) 

When three consecutive elements of a chain 
occupy the central site and two of the first 
shell sites, and the other elements of the 
same chain chain occupy none of the first 
shell sites, then the outer part of the local 
system weights the configuration with a factor 
Aone 

When the chain takes a configuration in 
which the chain passes through the central 
site and occupies /—1 sites out of the z—2 
first shell sites, then its weight is 2.,.; in 
particular, if the end element of the chain 
occupies the central site, the weight is ¢.2.,1. 

When the chain occupies one or two of 
the first shell sites (without occupying the 
central site), the weights are 2, and ¢,2 
respectively. 

For simplicity, the configurations in which 
the elements of the same chain occupy three 
or more of the first shell sites, without oc- 
cupying the central site, are not considered. 

When a solvent molecule occupies one of 
the z+1 sites, the weight is 4. 

The two functions A and B are defined 
by 


A = DAgKaa Se Asap 
B= gXap+ApXpp 


where 


Xin =eXD an } 


sn =exp( a } 


ope ( aon } 


The local partition function G is written 
as follows; 


G=EGat Gi, 2) 


where 


Coaldale ae )\@ietaat Ari 


a Zé 1ha, U 2 ee ‘)) CU Be INGE 
(2. 2) 


Guri=( 2 Viaye +é2-1ha,2-1A +2Ez2ha,z (2. 3) 


C= aZ(( : )) (Gieetag?)"BH, 2.4) 


The symbol (ay represents the 


number of ways to select r site pairs from 
among z—/—1 sites. 
Three cases; (I), (II), and (I); will be 


considered in the following. 


(I) Suppose 
Aav=0, t2=2 and ¢=0. 


This is the cass in which a chain never 
passes through the local system two or more 
times. 

Orr’s work‘ corresponds to this case. Ac- 
cording to his results, the average number of 
(a~8) contacts, 2X, is given by 


2Qa=nz—2n+2 (2. 6) 

B=1+441—6) et =1) CQ) 

—— (2.8) 
QaNat Ng 

2Qw =2éep— (Eaat+Epe)- (2.9) 


This is the well-known case in which the 
first order Bethe approximation leads to the 
same formula for X as that given by the 
quasi-chemical equilibrium. , 

Miller and Guggenheim‘ obtained the 
same result as Orr. 


(11) Next assume that 
jui*0, 121 and ¢=0. 
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The procedure of the first order Bethe 
approximation makes the probability of oc- 
cupation of the central site and the first shell 
site by a and § molecules respectively equal 
to the probability of occupation in the reverse 
manner. This equivalence is expressed, in a 
formal way, by the relation 


AG dG 
A =A 
Ody 


Oha 


Hereafter, in the equilidrinm state, the re- 
lative probability of the occurrence of the 


© 
O O 
(a,1) form Aa, (a,2) form da,» 


O 


(a,3) form ra,, (2,1) form ¢, hay, 


NI. 


Oh OM VO; 
oO 
a2) form e,Aa,, pha 
re el 3 
Cm® © © © 
O Oo 
wa Solvent molecule Ag 


Fig. 1—Various configurations and their weight 
factors in the local system. 


REVIEW OF THE ELECTRICAL COMMUNICATION LABORATORY 


configuration of a molecules in the (a,/) form 


(see Fig. 1) will be described by p(/). 


Ga;1 


2 
a (2. 11) 


=p). 


The probability function p(Z) depends on the 
state of the whole system, including the 
concentration and the temperature, and it 
may be considered to have a close relation 
with the equilibrium extension of the chain 
molecule. 

Furthermore, the ratio of the probability 
of occupation of the central site by the end 
elements to that of the other elements is 
represented by q(Z) as 


Rede pA! 
ie dart 
2 


From Eq. (2.10), the following relation is 
obtained by the use of Eqs. (2.2), (2.3), (2.11) 
and (2.12). 


=). (2512) 


v4 


ZpAghgXapB? = ( D 


\aasraeteg 


og) sii ie 
Dig PD {@—-l-1)+@—-Daq@}. 


(2. 13) 


The equilibrium properties of the whole 
system are also given by 


Dp. pao sar 
vere (2. 14) 
where v, and v; are the volume concen- 
trations of the a and the $ molecules re- 
spectively. 

In this paper, the following abbreviations 
are used: 


LO “Fe 
edt {2-Il-) +@—Dq} 
HCD) 
n 
_ nz—2n+2—9 


(2. 15) 
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e=2U-) 20,1 
i 


201,10 
m1 


an Oa 
1=2 


q)= 


> 


where Q;,, is the probability with which J/—1 
elements of the same chain, other than the 
two elements adjacent to the elements at the 
central site, also occupy the first shell sites, 
when the jth element of the chain passes 
through the central site. 

X is given by 


X 


} se OG/OXag 
NNg+Ng : 


2zG 


(216) 


From Eq. (2.16) the following equation is 
obtained 


= da NaNp 2 
= ( Pal 
Ga!NatnNp 148’ ake 
B2=14+4' (1—£’) (@"/FF—1) Cals) 
py GaNa (2. 19) 
Ja’ Na+Np 


zqx may be considered to be the effective 
number of possible contacting first neighbors 
of the chain molecule. Therefore, it is seen 
from Eq. (2.17) that in case J], the first order 
Bethe approximation and the quasi-chemical 
equilibrium are equivalent in some sense. But 
it is important that gq,’ should be determined 
so as to minimize the free energy of the 
whole system. gq,’ is never constant but de- 
pends on the concentration and the temper- 
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ature of the system. 


(II) In the most general case, the following 
relations hold: 
iim (St and 6-0. 


The relations shown by Eggs. (2.10), (2.11), 
(2.14) and (2.16) should also hold in this, 
case but instead of Eq. (2.12), the following 
equation is used: 


22th Z(( 2—l Vase dy Ar 
; =4()). 
6s: (ea) coma 
(2. 19)’ 
In addition to the above relations, the 
following relation is also necessary 
0G 0G 
ee Roe 
~ a0 p ae (2. 20) 


The above relation Eq. (2.20) show that a 
proportionality exists between the number of 
the chains which occupy only one of the 
first shell sites and that of the chains which 
occupy two of the first shell sites. Obviously 
R has a close relation with the equilibrium 
extension of the chain. If the chain spreads, 
R becomes larger; and if the chain contracts 
R becomes smaller. R, as well as e, should 
be determined so as to minimize the free 
energy of the whole system. 

With the help of Eqs. (2.10), (2.11), (2.12’), 
and (2.14); Eq. (2.16) gives the expression 
for X as, 


XX” 
Zh 
a! Ng ——* n\n 
aG <a ae Sie een as ae 
Ja!" Na+Qp’ Np _— 4Arsnp/Z ae, Aga” NaNp (e2v/KT 1) 
Qa" Nat Qa’ Ms (a"NatQp'Ns)” 
2rs ngrere /KT 
x” —— é —e = —— = 2 SS 7 = = = 
a" Nats’ Ns eS 2ng __(e2wiK? 1) ren Aga" NaNg —(e2wiKT _}) 
Qa! Nat+Qe Ng Ga” Nat+Qp Np)” 
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206 ow /KT 
OU aaa ak a aes 


where, re and 7; represent the numbers of 
chains which occupy two of the first shell 


sites when a and § molecules occupy the 
central site respectively. 


1 


| pees 


- rel 
ai aa 


=( Zz i Gea )\rprtradt art 
Ca l= 2 if 


ayl 


= i ( : )) Ace es 


i= CG 


To obtain r, and 7, we must make use 
of (2.20). Unfortunately, this prpcedure is 
so complicated that it seems to be difficult to 
obtain the results in a closed form. 

For example, for a two dimensional hex- 
agonal lattice, if the chain is sufficiently long, 
so that 


<1 
n 
we obtain 
Ay 
ey ic {jee 
n n , 
T2—= RyQ Yt R+2 —VaNsreres 


This problem will not be considered further 
for this case, but it is noteworthy that in 
case (III), the equivalence of the first order 
Bethe approximation and the quasi-chemical 
equilibrium are no longer correct. 

In the following, we confine ourselves to 
case (JI), because the formulas obtained in 
this case are relatively simple and yet seem 
to retain the desirable characteristics. 


3. Thermodynamical Quantities 


The partition function of the whole system, 
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Q, is written as follows: 


Q=>2(p). B.D 
pg 


The partition function for fixed p, 2(0), can 
be written as 


/ —2X(p)w | 

2 =2 = ne ay n i ee x 
(p) vexp| KT Or . py P, X) 
Gr 

Ea) —l Zz 

Qo =CaMa«Cs"2EXP Gea t) 
(353) 
fine ab Xp) a I (3. 4) 

0 


where, ¢, and ¢; are the kinetic energy terms, 
Yw(p, X) is the number of configurations for 
x 


constant p, zX(p) is the statistical mean of 
zX for the constant e, and 


aX (p= = eTe ee 
Ow 


We can generally follow Orr’s method‘ 
in deriving the terms X and Yo(po,X). Orr 
used z(z—1)"-2 for the number of configu- 
rations of the isolated chain. This factor is 
replaced by /f(e@) in our case. /f(o) is the 
total number of configurations which the 
isolated chain of constant po takes. 

The derivation of f(e) will be deferred 
until the next paragraph. 

The free energy of the whole system, F= 
—kT log 2 can be approximated well by the 
maximum (0) in the summation of Eq. (3.1). 


F=—kT log Q(p). 


The following formula of F is obtained by 
the formal extension of the Orr’s formula to 
our case. 


ee kT p= Jee" 
W ae —0 i ae ag eal See 
aD (a Oe 
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/ _9¢E/ 
Bo log B’ +1—2€ 


Gs) 1--) 


seg ee ya” 
n-- | nde n— Qa’ 


2 ue n n-i 
Pall AR Rl ae a eS 
bia f(p)z(z—1) ea 


n ; n 


—0,)+0, log ve tnd —v.) log d-»)}, 
55) 


The free energy F must satisfy the re- 
lation 


or 


SRT. 10 ae rae 


f n—1 = n- ) n 
ae pele 
nz Nn— Qa N— Qa 


ee asian! bu! a 
0 Zz 


2 
ee i 7) tog (<7 -2) J=2 
NZ N—Qa N— Qe 


(3.6) 


cam RT 


op is determined from the above equation. 
The mixing energy AH is shown by 


yee ie a. i : 
da Na tne “148 


Garo) 
The mixing entropy AS is derived from 
TAS=AH— AF. (3. 8) 


The chemical potentials of the a and § com- 
ponents are 


1 1 a ai aE) do 
Hie ee PT Cie Jig, RT X80 Jnany 0 


=e 
eet ae 


Ti pe ZOR VS Wie eee 
fap faa 5 —lo 
<n i SMoRGs el, gate Picea 
_ M1) qe" n 
=| - ib GaP IES = au 
on (5 nf(o) + pa log ee 
SO aM reals 
N—Ga! : | 2—qa! | tog re 
3. 9) 


a Br 1-26" 
er =~ opptety 8 Gas wy 


E(t st 
(3. 10) 


The critical mixing point is determined by 
the following simultaneous equations. 


Ofte Oe. 
OVe Si Ov—2 ay 


The resultant formulas which give the 
critical point are very complicated, but if the 
chain is sufficiently long, it is found that the 
critical concentration is proportional to 1/7, 
and o is approximated well by linear form 
of v, at that concentration. 


(Va) = Pot pide. 


The approximate solutions of the critical 
point are 


/ 


yogi 2 
base Tn V (2—2—po/m)(2—1—2 p/m) +1. 5p1/nz 
Ga) 
2w 1 2 1 
exp—__=1 


Els, . 2 -2aph | eet ey 


ZW 2p _ Pie 
a) 7-2—pjn | Ga2=o/nye xn 
(Gy IW) 
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4. The Derivation of f(p) 


Perhaps the most orthodox method to ob- 
tain f(e) is a matrix method. The method 
presented is a matrix method modified slightly 
to obtain the partition function of an isolated 
chain molecule. By a step-by-step process on 
the lattice space the configurations of a chain 
are constructed from one of its ends. When 
the rth step is finished, an allowed configu- 
ration of m segments from the (r—m)th 
to the step is called o(7). The interactions 
between the segments more than m steps 
apart are neglected. The number of all 
possible configurations of m segments is ». 
These configurations are designated by ¢;,(7), 
O(N), y.(r). When one step is added, any 
configuration of m segments from the (r—m 
+1)th to the (r+1)th correspond to one of 
these y,;; these are designated by 9,(7+1), 
@2(7-+1), v0 g.(r+1). The factor 7 is weight- 
ed to each process, if the segment of the 
(r+1)th step has » nearest neighbor contacts 
with some segments from (r+1)th to (r—m 
+1)th. The weight is 0 when overlapping 
occurs, and 1 when there is no interaction. 

In matrix form, this is described as 


AaAr+)D\ =) (aM (4.1) 
g2(r+1) ¢2(1) 
g.(r+1) o(1) 


7ij=0, 1, 7, +--+ ne}, 


The largest eigenvalue of (7:;) is an(y), 
and then, for sufficiently large 7, 


p(n) =” om). (4. 2) 


On the assumption y=e?, 7°” is the 
partition function of an isolated chain. 

When 2”,,”” is expanded into a power series 
of 7, f(o) is the coefficient of yl. If f(p)yrl 
is sufficiently large for some value of 7 and 


n>m, 


log f(o) =n: log in-(5) logy, (4.3) 
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We consider that this p is equal to its 
mean value 


Py F108 dm 4.4 
mi rp (4. 4) 


Differentiating Eq. (4.3) with respect to 9, 
we obtain 


d It 
ce log =. log 7. Ga) 


The above relation must be substituted into 
Eq. (36). Particularly, if the concentration 
of the chain molecules is sufficiently dilute, 
Eq. (3.6) is reduced to 


w d 
2 tog f=6 
eT! dp eet 
or 
n= ew kr 


This is a satisfactory result, because in this 
case 4” must be the partition function of 
the chain. 


5. Numerical Calculations 
5.1. The Mixing Energy is Expressed as 
AH=zXw. GD) 


In case (] ) (Orr, Miller, and Guggenheim), 
X is given by Eq. (2.5). According to this 
equation, the curve of AH versus €(1—6&) is 
expected to be a single curve. Contrary to 
this expectation, however, the measured curve 
of the mixing energy of polymethylmetha- 
crylate and trichloroethylene against &(1—6) 
is found to be not a single but a closed 
curve.°” (See Fig. 2.) 

This is explained by Eq. (2. 17) in case (IL). 
Because gq,’ is not constant: but changes in 
value as the concentration increases, the curve 
AH versus &(1—€) is never a single curve, 
it is a closed curve. 


5.2. Critical mixing Point 


In case (II), the critical point is given by 
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io 


Py ee) 
@ 10 (> 
of 
& 
S 
a Q is given by 
= 
a 5 
<<) n 
QD) =2", 


where 4 is the maximum eigenvalue of the 
following matrix: 


0 0.1 0.2 0.3 (4 a) 
\ ds) 
Fig. 2—Value of AH/RT(qana+nzs) plotted 
versus €, 6:(z=6) for solutions of Wee 14 
_ polymethylmethacrylate in trichloro- ie , 
ethylene; X: experimental points. 
Data from Staverman and Dekking.“ ne=1. 74 
Po_—0, 646 
Eq. (3.11) and Eq. (3.42)... These formulas ano 
are reduced to those of Guggenheim when ie 
Po=01=0. So, in case (II), a larger value 
: : slOWl = 
of vz, can be predicted than in case ([). ee are 


Oo and p; in case (II) are given by 
whereas Guggenheim’s formula gives 


d 
po=2r( 4 oe a) 
me ee pee 
COS Cie Vn ’ 


dlog Q ad’ log Q 
=2 {( be ( ) | 
a | dy eee dy? va=0 and Flory’s theory gives 
Y ae 2(42—2—o/”) Q—-1) | Uxer—l/ 9 
nz Z 


(5.3) 5.2.2. Three Dimensional Cubic Lattice 


The “fundamental vectors” are 


The following are numerical calculations 
for simplified cases. 


5.2.1. Two Dimensional Square Lattice 


In deriving Q(y), the “fundamental configu- 
rations” or “vectors” are 
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A is maximum eigenvalue of the following 
matrix: 


uff), Su brea een f 
Te eels wad 
ib at AEP 
jan 7 0 
g— 3+ V9+129 
2 
iP _—(), 391 
n 
-P1 — _0, 032 
nN 
fh oe 
aye Vn’ 


whereas Guggenheim’s formula gives 


Pens 
aCe — Vn , 


and Flory’s theory gives 


fen 
OTN OY le a Vn 


The experimental results‘? show that 
Seesaw Vn. 


Therefore, our formulae Egs. (3.10) and 
(3.11) predict a better value for v,,. that the 
previous formulas. But the agreement is not 
yet satisfactory. 


7. Discussion 


It is already known that the binodial curve 
of the phase diagram of the polymer~solvent 
system predicted by Flory’s theory is too 
narrow to explain the experimental results. 
The author has not yet calculated such a 
binodial in the present paper; however, he 
considers that a broader binodial will result 
from his theory than from Flory’s. 

This is reasonable because, in our theory, 
the concentration dependency of the extension 
of the molecule is considered. 

It has been found that our theory can 
explain the experimental facts better than the 
previous theories, but the numerical calcu- 
lations given in the preceding paragraph 
show that there still remain some differences 
between our theory and experiment. The 
reason is partly because of the fact that the 
number of the “fundamental configurations” 
or “vectors” g, used in calculating Q(y) is 
too small. In the calculation for the cubic 
lattice, only ¢, has an intrachain interaction. 
This means that in the summation era 


1) 219i, the terms for 7>3 are all neglected. 
Jj 


If all of these terms are included, o becomes 
larger, which makes v,,. larger. 

Even if the factor eo is correctly taken into 
the calculation, some unsatisfactory points 
will still remain. For example, the curve of 
AH versus &(1—€&) may not be expected to 
be such a large loop as the experiments 
show. However, if the theory can be con- 
structed from case (III), further improvements 
can be expected. Unfortunately this seems 
to be very complicated. 

In addition to the above discussion, the 
concentration dependency of p will yield 
valuable information concerning the concent- 
ration dependency of the equilibrium ex- 
tension of the molecular chain. 


7. Conclusion 


By using the first order Bethe approxi- 
mation, the author tried to include the effects 
of the extension of chain molecules in the 
theory of the solution of chain molecules, 
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and calculated the various thermodynamical 
quantities. It was found that the effects of 
the extension of chain molecules are important 
to the problems of the mixing energy and 
the concentration at the critical mixing point. 
Our theory applied to those problems gave 
better agreement with experimental results 
than did Flory’s theory. 

Further, it was shown that the well-known 
theorem of the coincidence of the result of 
the first order Bethe approximation with that 
of the quasi-chemical method does not always 
hold in the solution of long-chain molecules. 

The author considers that further efforts 
to include the complete effects of the extension 
of chain molecules in the theory will yield 
more satisfactory results. 


(1) 
(2) 


(3) 


(4) 


)) 


(6) 
(7) 
(8) 
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A Channel Selector Applying Error- 


Correcting Codes’ 


Shoji YOSHIDAt and Keisuke TOMARUt 


This paper describes a channel selector, which makes use of error-correcting codes. Para- 


metrons are used as the source of selection signals. 


The operation of the selector is stable, 


and its time for selection is short. It can be constructed simply and economically. Causes 
of variations of the selecting ratio for a partially equipped selector are clarified, and 


some countermeasures are suggested. 


1. Introduction 


There are many types of channel selecting 
circuits, such as diode matrices and a tree 
circuits, that are used in exchanges and com- 
puters. The function of these circuits is to 
select one specific circuit or element out of 
many. 

New type of selecting circuit is presented. 
Nature of error-correcting codes is applied to 
the selective principle of the selecting circuit. 
This method enables a simple and economical 
construction of the selecting circuit. Moreover, 
the time needed for its selective action is 
usually shorter than that of a tree circuit. 
Parametron signals are used as input signals 
in this selecting circuit. Oscillation phases 
of parametrons (0-phase and z-phase) are 
made to correspond to the “1” and the “0” 
respectively of binary selecting codes. Small 
ferrite troidal cores with a linear character- 
istic are used as circuit elements. This kind 
of selecting circuit can be widely used. 


2. Selective Principle 


For any two codes with the same code- 


MS in Japanese received by the Electrical Communi- 
cation Laboratory, June 8, 1960. Originally published 
in the Kenkyi Zituyoka Hokoku (Electrical Com- 
munication Laboratory Technical Journal), N.T.T., 
Vol. 9, No. 6, pp. 585-620, 1960. 

Switching Research Section. 


—+ 


length 2 (the number of bits in the code) 
and with mutual destance d, the selecting 
ratio R; is defined as follows: 


n 


de @ 


s 


Now let us consider the following two 
codes: 


X= Cp Xo, crrrttees ; Xe) 
(2) 
Y= y, Yo, ttttesees EE 
d is given by the following equation: 
d=dis(X Y)= 3 (wy)? (3) 


For example, if the codes X and Y are: 


A=l 0 0 O30" 2 
4 
¥=0) 170. 07 0ar 


d is 2. (n=6 in this case.) 

If all the code-lengths are the same and 
the mutual distance between any two codes 
is not smaller than dy, the following equation 
can be obtained: 


n 
sa : 
n—2dp ©) 
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By making R; larger, higher selection stabili- 
ty and higher selection reliability can be 
expected. 


Table 1 


SELECTING CODES 
Nelo" 7=6, B=3) 


a . 
i 1 ao eG) TG 
0 (= 0 0 0. 0 «0 
1 (PrOe sre = OA 
2 (eer we. 0 th 
3 ie Mie ee OTS nO 30 
4 Gees ee sO 
5 PO Ona St 
6 Ty TO ee OS ee es | 
7 ite Ce et a 
8 ets ey a ee aw 
9 eee! ore Nal 
10 He ieee A yt 
11 Vette Olea At 6 
12 (seein 1-0; .. 0 
13 os, Me ee eh a ee 
14 Oeste i Or © I 
15 ‘OS eee ee ae 


Table 1 indicates the selecting codes for 
n=6, R;s=3 and N=16. (N is the number 
of the channels to be selected.) Correspond- 
ingly to this table, six input windings are 
wound on each core, as illustrated in Fig. 1. 
Each core has an output winding with a 
tuned circuit. The output voltage vs. input 
ampere-turn characteristic of each coil is 
linear as shown in Fig. 2. 

When six input currents are applied simul- 
taneously to the windings in the directions 
indicated in Fig. 1, magnetic flux is induced 
in each core in the direction shown by the 
arrow. If the value of the magnetic flux 


Output 


Input windings Cores Circuits 


(ee Sars 


y induced flux 


i| of of of of i| | trout current 


Fig. 1—Construction of selecting circuit. 


lin 
[vou 
Vout 


Vo 


ATin 


Fig. 2—Characteristics of the selecting-circuit core. 


induced by each input current is equal to J, 
a total magnetic flux of 6J is induced in 
core #1 in the same direction and an output 
voltage Vo is induced in the output winding. 
As the value of the magnetic flux in core 
#0 or in core $15 is only 4J—2J=2J, the 
consequent output voltage is only V)/3, which 
is smaller than that of core #1. In this case, 
the selecting output voltage appears across 
the output winding of the core #1, for the 
code pattern of input current and that of 
input windings of core #1 coincide with each 
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other. If the load circuit operates for the 
voltage Vy while it does not operate for 
V,/3, the selective operation is performed 
satisfactorily. 

Since the oscillation phase of a parametron 
is either in phase or opposite phase to the 
stanpard phase, patametrons can be used as 
the input power source of selecting signals. 
For reference, the relation between the code- 
length (7) and the number of the channels 
to be selected (N) is shown below, when the 
selection ratio is 3. 


Ah dee ee) 
Ni4 16 64 256 


3. Design of the Selecting Circuit 


3.1. Characteristics of the Ferrite 
Cores to be Used 


In the design of the selecting circuit, the 
values for the following items are given pre- 
viously: 

1) Output power. (Output voltage, load 
impedance, etc.) 

2) The minimum selecting ratio. 

3) The number of channels to be selected. 

In order to obtain the required output, it 
is necessary first to select a code which 
satisfies the above conditions, to determine 
the input current value, to determine the 
turns ratio of windings, etc. Special attention 
should be paid to the selection of the cores 
to be used, because saturation of the core 
will reduce the selecting ratio. Input vs. 
output characteristics of some ferrite cores 
are illustrated in Fig. 3. 

In the subsequent experiments, L,;D,, double 
cores, made of Cu-Zn ferrite for high fre- 
quency use, are used to avoid the reduction 
of the selecting ratio caused by saturation, 
and to obtain a large output voltage. 


3.2. Selection of the Core 


Variations in the output voltage between 
cores cause a reduction of the selecting ratio. 
Therefore, variation must be kept to a mini- 
mum. 
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R: 10KQ | 


Vout( Vp ra o) 


lin (AT, )— py) 


Fig. 3—Input vs. output characteristics of 
ferrite troidal cores. 


In the secondary tuned circuit which is 
illustrated in Fig.3, the variation of the 
output voltage depends mainly on the _ vari- 
ation of Q of the core used. Therefore, QO 
was measured to provide a criterion for 
selecting the cores. The correlation factor 
between Q and output voltage was 0.85. An 
allowable output voltage variation of +5% 
was assumed, and cores which had a Q= 
26~36 at 1 Mc/s were selected. This mea- 
surement was carried out after winding a 
ten-turn coil around the L,D,, double cores. 


3.3. Design of the Input Coupling 
Circuit 


As illustrated in Fig. 1, the input windings 
of the cores are connected in series. If the 
output circuits are complete, the mutual 
couplings between the input windings are 
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negligible, because of the nature of the code. 
oo the input impedance is given 
Wi 


Li nR. (6) 
Core: L, D,«(double ) gat SOT 
Ro=50K IN solid lines Ty ( R, 
30 -R,=30K Qin dotted lines M 
| 
20 + MS C — 


R, X(Q/1Me/s ) 
iS 
\ 
ee 


[0 80 30 100 MPU oe~ 120 190 
y 1 Or) 
10 ya 
| 
| 
20 : 


Fig. 4—Input impedance of the selecting circuit. 


R is load resistance seen from one input 
winding of a core. The value for R is given 
in Fig. 4 as a function of the tuning capa- 
citor. The power (P) consumed in the load 
is given by, 


P= le @) 


I is the input current required to obtain the 
required output power. Input signals must 
satisfy all these conditions. And then in 
practice it is not desired that the input current 
is varied by the load effect, for the variation 
reduces the selecting ratio. Under a partially 
equipped condition, the input windings couple 
among themselves and reduce the selecting 
ratio. For these reasons, it is desirable for 
each input circuit to be a constant current 
source. On the other hand, the operation of 
parametrons becomes unstable with a increase 
in the value of the load impedance. This 
instability makes the use of parametrons as 
input signal sources undesirable. 


Parametron Grounded — base 


unit A class amplifier Selecting 


circuit 


Low imp —> 


<— High imp 


Fig. 5—Transistor coupling circuit. (Impedance 
converter.) 


Because of the above reasons, an impedance 
converter is required which has a low input 
impedance to ensure stable parametron oper- 
ation and a high output impedance to operate 
as a constant current sourse. Such an im- 
pedance-converting input-coupling circuit is 
shown in Fig.5. A common-base amplifier 
using a Germanium transistor with an alpha 
cutoff frequency of 10Mc/s is used. The 
input impedance of this circuit is approxi- 
mately 25Q and the output impedance is 
approximately 3kQ at 1Mc/s. Therefore, 
the constant-current-source requirement is 
satisfied. Besides, the operation of parametrons 
is not disturbed by load effects. The main 
reasons that the transistor is employed in the 
common base configuration are as follows: 

1) The largest output impedance may be 
obtained. 

2) The influence of the variation in the 
current amplification factor of the transistor 
is less than that in other configurations. 


4, Experimental Results 


Experimental results are tabulated in Table 
2 and in Table 3. The cores used in these 
experiments were selected according to the 
principle explained in the preceding section. 
The experimental result shown in Table 2 
corresponds to the error-correcting selecting 
codes presented in Table 1. The operation 
of the selecting circuit was examined in both 
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fully and partially loaded conditions. The 
selected ontput voltage was approximately 13 
volts. 

In the experiment whose results are given 
in Table 3, another code (N=256, n=18, 
R;=3 or 9) was used, and 48 loads were 
partially loaded. In this experiment the 
selected output voltage was approximately 18 
volts across a 50kQ load resistor. 

In Table 2 and Table 3, the term “direct 
coupled” means that input windings of the 
selecting circuit were connected directly to 
the output terminals of the parametron units, 


Table 2 
MEASURED SELECTING RATIO 
CNEN6weRs—3)) 
No. of equip- Input con- | Selecting ratio 
ped circuits nection : 
| mean max min. 
Direct 
16 (full) cacpled Di 2.8 2.6 
Direct ee 5 if 
12 (part) boavled 2.8 yo) Za 
TR 
12 (part) coupled peo 3.0 Oa 
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and the term “TR coupled’ signifies that 
common-base transistor amplifiers were insert- 
ed as illustrated in Fig. 5. 

In Table 3, the symbols “+” and 
are used as follows: When a code pattern of 
the input current is compared with the pattern 
of any input winding, and when correspond- 
ing digits are the same, this condition is 
represented by “—”; whereas when corres- 
ponding digits are opposite, this condition is 
represented by “+.” After all the digits are 
compared, the number of “+” and the 
number of “—” are added, and when the 
former is larger than the latter, this condition 
is indicated by. “+” in Table 3, and the 
opposite by “—.” 

For example, in Fig. 1, the code patten of 
the input windings on core $l and the pat- 
tern of the input currents are the same. But 
at core #0, they are not all the same: 4 same 
digits and 2 opposite digits are seen. There- 
fore, this situation is represented by “—.” 
In the same manner, the symbol “+” is 
given to core #15. The difference between 
“1+” and “—” will be explained in the 
following section. 

Oscillation phases “0” and of the 
parametron are made to correspond to “1” 
and “0” respectively of the binary number. 


“ce bed 


eo 9» 


T 


Table 3 


MEASURED SELECTING RATIO N=256 (Number of equipped circuits: 48) R;=3 or 9 


5555 EE eee 


oe . 
Vie Selecting ratio 3 9 
I Hs ; et 

nput connection ye Oh na max. min mean max. min. 

: ja : 
mean Bh bul Dt ZS 50 4.4 
Direct coupled th 3.5 5.4 29 | 24 50 9F3 
- 2.2 nee Pa oy, 5) 7.8 4.4 
mean 2.8 3h 25) Tay | 10. 2 4.5 
TR coupled ate 2.9 3S} 2a 8.7 10.2 50) 

=). ee eee - Nees (0 
= 23 rae 8) 213% | 5.6 6. 2 4.5 
| =| | 
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The variation of the input current due to 
the variation of the oscillation voltage of the 


parametrons and to a of the transistors was 
smaller than +5%. 


5. Variation of the Selecting Ratio 
5.1. Causes of Variation 


The following terms are considered as the 
causes of variations: 

1) Variation between cores. 

2) Variation in the input current. 

3) Phase shift between input cutrents. 

4) Asymmetry and distortion of the input 
current wave form. 

5) Change of the input current due to the 
effect of coupling between input windings. 

The effects of 1) and 2) have deen de- 
scribed. Term 3) is negligible when para- 
metrons mounted on the same wiring board 
are used. Term 4) is also negligible when 
parametrons are used as input signal sources. 
Therefore, only the effect of 5) is thought to 
be large. The coupling between input wind- 
ings takes place mainly through the core as 
an intermediary. The effect of the winding 
stray capacitance may be neglected. 

When all the cores and the load circuits 
are mounted, the coupled current flowing in 
each input winding is completely cancelled 
by the characteristic of the code. (It is as- 
sumed that coupling factors on every core 
are the same.) However, when the selecting 
cores are partially equipped, the values of 
the input currents will be changed because 
of the incomplete cancellation of the coupling 
currents. Therefore the variation in the 
input currents may become large. 

In order to prove the above assumptions 
and to explain the causes of the variation of 
the selecting ratio indicated in Table 2 and 
in Table 3, some basic experiments, as de- 
scribed below, were carried out. 


5.2. Basic Experiments and Consider- 
ations 


A selecting circuit, illustrated in Fig. 6, 
was constructed using the code indicated in 


Input windings Core Output circuits 


1 2 3 
9 357 
0) 0 0 
: (0) 50KQ. 
1T 17 


oh 
b 


Fig. 6—Configuration of selecting circuits. 


Table 4 


SELECGLING CODES CV=4, n=3, Rs=3) 


—— C | 
SS | 1 Zi 3 
0) 0) 0 0 
1 1 0) 0 
2 0 il 0 
3 1 iL 0 


Table. 4>(N=4, n=3, & Rs—3.) The coupk 
ing between the input windings 1 and 2 
through core #0 is indicated by a2, and the 
other coupling constants are indicated in the 
same manner. Input currents without coupl- 
ing are presented as J,, J2, and J;, and those 
with coupling are presented as J,’, J;', and 
I,'.. By taking into consideration only the 
coupling current of the primary term, and 
neglecting those of other terms, the following 
equation is obtained: 
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ie | 1 DiGi Di@13% || 2s | 

| 2! | Die@e1i il Yaoi | I (8) 

iT; | | Masr = as2i 1 || Zs 
DAG. e. 


Now the following conditions are assumed: 
1) W=|h=|i=L, (9) 
2) The absolute value of @ is constant. 
3) The sign for a is “—” when any pair 

of the windings have the same directions, 

and “++” when they have oposite directions. 
4) The output voltage is proportional to 
the total sum of the input currents applied 
to that core. 
These assumptions are appropriate for prac- 
tical purposes. 


(A) Fully loaded condition (See Fig. 6) 


Equation (10) is obtained from Eq. (8). 


| Li’ | [1 0 O||Al [AI 
ES i) Wel ava (10) 
Ca OrenO meat ep LT| 


The values of the output voltage of each 


Tabie 5 


CALCULATED & MEASURED OUTPUT 
VOLTAGE (VOLT) 
(Under the fully equipped condition) 


Ine 1 Deel ges: 
_Core number ~ | | 


31 1 I I 
0 
8.0 he Hh 2.8 De Ut 
I 3L I 
wt i 
Dill oe 2.8 2.8 
I If 31 
. I 
2.8 Pi) 8.0 2.8 
I Il if 
3 or 
2.8 208 Zell G8 
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core are indicated in Table 5. Input currents 
are shown above the oblique lines. Their 
values are proportional to the values of the 
output voltage, as explained above. Measured 
values of output voltages are shown below 
the oblique lines. In this experiment as well 
as in the following ones in this section, input 
windings of the selecting circuit were coupled 
directly to output terminals of the parame- 
trons. The variation in output voltage shown 
in Table 5 was mainly because of variation 
between cores. 


(B) The case in which core #0 
is omitted 


Equation (11) is obtained from Eq. (8). 


Te 1 acre cma 
hRY\=|a 1 a@ | ld: e5 
Ihe ee se NIE 


Two digits of the binary code, “0” and 
“1,” were made to correspond to “+” and 
“_”” respectively. For example, if the input 
current pattern is (1 0 0), 


=e and JE — FN b 


Therefore, from Eq. (11), 
I,’ =—f+2el,_ and f,’=],’=f. 


Since the input winding pattern of core #1 
is (1 0 0), the total input current given to 
this core is 


Ti! +1! + 2,’ =I—2al+I+1=3I—2al. 


On the other hand, as the input winding 


pattern of core #2 is (0 1 0), the total input 
current is 


Rly See PS 
=—(I-2al). 


In Table 6, only the absolute value is 
taken into consideration as the output voltage, 


and both calculated and measured values are 
tabulated. 
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Table 6 
CALCULATED & MEASURED OUTPUT VOLTAGE (VOLT) 
(Under a partially equipped condition I) 
ee See 2 eee 
Input code ] 
: 0 | 
Core number | : | g : 
il I+2aI 3.0 3I—2aIl 8.0 | I-—2al 2.6 I-—2aI 2.6 
2 I+2aI 3.0 I-—2aI 2.6 3I—2aI 8.0 I-2aI 2.6 
3 | T+2eF—3.0 jh, Cae ee ye On 31-2aT 8,0 
| / 
Table 7 
CALCULATED & MEASURED OUTPUT VOLTAGE (VOLT) 
(Under a partially equipped condition II) 
Input code | 
0 1 2 3 
Core number | 
2 ie 258 Ik Baits 3l—4aI 8.0 I—4aI 2.5 
3 JE BESS IE ss I—4aI 2.5 3l—4aI 8.0 
(C) The case in which both cores #0 
and #1 are omitted ; | 
ES | CG jest OV5=p | 
a ae 
Both calculated and measured values of the - | _Aboazsy ia 
output voltage for this condition are tabulated balks | [_iMog 5 ao 
in Table 7. 3 aes = er 
: ~ é yes oo | 
Since the measurements for cases (A), (B), S si / ia enn. 
and (C) was done independently, it is mean- 3 | | eat “| 8:28(at_IMe/s) |. 
ingless to compare the-selected output voltages Ss jt a Mest a fo 
of each case with one another. pa - wile —\\ Coupling factor = Wet 0.2 5 
Se 
(D) Coupling factor a ae are ie a et SS Se 3 
—— | | j | | | = J 
‘Gen 2 40 60 80 Too 120 
In the above experiments, the coupling C(pF) 


between any two of the windings is consider- 
ed a result of the effect of the output wind- 
ing. In other words, if some of the output 
windings are opened, the same results are 
obtained as those when the corresponding 
cores were not equipped from the beginning. 

Fig. 7 illustrates the measured values of 
the coupling which occurred between two 
windings due to the effect of a tertiary 


Fig. 7—Change of the coupling constant due 
to the effect of a tertiary winding 
with a tuning capacitance. 


winding with a tuned circuit. The core was a 
L,D,, double core (Q=28 at 1 Mc/s). When 
the tertiary winding was opened, namely 
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when C=0, the coupling constant between 
two coils with a one turn winding each was 
approximately 0.006 (44 dB). This coupling 
constant is given as the ratio of Vout to Vin 
in Fig. 7. With an increase of tuning capaci- 
tance, however, the coupling constant changed, 
as is illustrated by curve 1. Curve 2 repre- 
sents the induced voltage in the tertiary 
winding. In an actual circuit, a resistance is 
added in parallel with the tuning capacitance. 
This causes a change in the coupling constant. 
Besides, the coupling constant is changed by 
detuning or the variation in the element. 

From Table 6 and Table 7, a was calcu- 
lated to be approximately 0.025, under the 
assumption that the output voltage is pro- 
portional to the input current. In the above 
experiments, the tuned circuit voltage was 
observed with an oscilloscope while the tuned 
circuit was being adjusted. When the oscil- 
loscope probe is disconnected, the resonant 
frequency shifts because of the removal of 
the approximately 12 pF of prove capacitance. 
This point is indicated in Fig. 7 by the small 
circle P. At this point Vout/Vin is nearly 
equal to 0.025. This value corresponds to 
the value obtained above by calculation. 


5.3. Other Experiments 


Finally, the calculated results of input cur- 
rent and measured values of output current 
under the conditions indicated in Table 3 are 
tabulated in Table 8: N=256, 48 cores, and 
output circuits were partially equipped. In 
this experiment, transistor coupling circuits 
were used. 


6. Some Hints to Improve the Selecting 
Ratio 


1) Use cores which have a large saturation 
flux density and good linearity. (In general 
these requirements are opposite to those for 
miniaturization.) 

2) Make the input signal source impe- 
dance high, avoiding variation in the input 
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Table 8 


CALCULATED & MEASURED OUTPUT VOLTAGE 
(N=256, Number of equipped circuits: 48) 


$e 


Input | Output vol- | 
Core eral tage (volt) | Input current 
Selected alts) | 18.5 18J—352al 
red 6.6 61 352al 
—6 8.0 61+ 224al 
Non- ¥ 4 
selected | 
+2 ZnO ) 2I—288al 
—2 3.0 21+ 288aI 


current due to the coupling. 

3) Clip the output voltage at a proper 
level. For example if the selected output 
voltage is 20 volts and the non-selected out- 
put voltage is 10 volts, the selecting ratio is 
2; whereas by clipping each output voltage 
at 5 volts, the selecting ratio becomes 3. 

4) Reduce variations between circuit ele- 
ments and between input currents. 


7. Conclusion 


By applyjng error-correcting codes to a 
selecting circuit in which small ferrite cores 
with linear characteristics were used, a good 
selecting circuit with high speed operation 
was obtained. 

The causes which lower the selecting ratio 
and make the selection unstable were explain- 
ed, and some methods of improving the ratio 
were suggested. 

The partially equipped (48 out of 256) 
selecting circuit has been employed in an 
experimental all-electronic telephone exchange, 


and has been operating stably for more than 
one year. 
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Reflections and Mode Conversions in Imperfect 
Junctions of TE,, Wave Transmission Line’ 


Ken-ichi NODA{ 


Imperfect junctions cause attenuation increase and delay distortion in the TE): wave 
transmission line. It is known, that the magnitudes of reflections at a slight tilt and at a 
offset do not have a first order term of tilt angle or distance between iwo waveguide 
axes respectively. Reflections at a tilt and at a offset have been calculated from the 


continuity condition of the electric and magnetic fields at the junctions. 


Mode conversion 


at the junction of true circular with elliptically distorted waveguides has also been calcu- 
lated by a similar method. The results of these calculations provide fundamental data to 
determine the dimension tolerances and installation tolerances for this kind of waveguide. 


Introduction 


Imperfect junctions in a TE); wave trans- 
mission line cause reflection and mode con- 
version. Imperfect junctions are classified into 
the following types: tilt, offset, diameter dis- 
continuity, elliptically deformed discontinuity, 
etc. Several of these types may occur simul- 
taneously in a practical case. First order 
solutions were obtained for the cases of the 
slight tilt-and the slight offset in terms of 
the tilt angle and the offset between the two 
waveguide axes respectively.‘ It was also 
shown that in both cases the reflection co- 
efficient does not have a first order term. 

In the case of diameter discontinuity, re- 
flection and mode conversion were calculated 
in the first order term of the diameter di- 
fference.°” 

Mode conversions due to imperfect junctions 
cause serious degradation of the characteristics 
of transmission lines. However, as the re- 
flection coefficient does not have a first order 
term except in the case of diameter discon- 


* MS in Japanese received by the Electrical Communi- 
cation Laboratory, on 19, July, 1960. Originally 
published in the Kenkyii Zituyéka Hokoku (Electrical 
Communication Technical Journal), N.T.T., Vol. 10, 
No.1, pp. 85-94, 1961. 

Hyper-frequency Research Section. 
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tinuity, its value is usually quite small. 
However, when a helix waveguide is com- 
pletely developed and used, the mode convers- 
ion may become little serious and the reflect- 
ion remains as a important problem to deter- 
mine the characteristics of this kind of the 
waveguide. 

In this paper; first, the reflection at a 
slight tilt, and second, the reflection at a 
slight offset are calculated from the continuity 
condition at the junction plane. The results 
are formulated in the form of a summation 
of the series of functions with regard to the 
eigenvalues of the modes converted at these 
junctions. 

The mode conversion at an elliptically de- 
formed discontinuity is an important value 
useful for determining the dimensional to- 
lerances of the waveguide. In a manner 
similar to that used for calculating reflections 
due to slight tilt and offset, the mode con- 
version can be obtained from the continuity 
condition of the electric and magnetic fields. 
The wave function in the elliptic waveguide 
has been obtained by the superposition of 
those of the circular cylindrical wave.‘ 

If the dissipation loss of the waveguide is 
disregarded in the analysis, the mode con- 
version factor increases infinitely when the 
frequency approaches the cutoff frequency of 
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some converted mode. Detailed analysis must 
be made in such a case to know the exact 
conversion factor. Taking into consideration 
the attenuation of the waveguide which has 
thus far been disregarded in the calculation 
of the mode conversion, it is found that the 
mode conversion at the cutoff frequency of 
a converted mode does not attain a very 
large value. 


1. Reflection due to a Slight Tilt 


Both ends of a waveguide are usually 
finished by attaching a flange or cutting a 
thread to connect to another waveguide. The 
end plane of the waveguide must be finished 
exactly perpendicular to the guide axis. A 
waveguide having a perfectly finished end 
plane and a waveguide having an obliquely 
finished end plane are assumed to be connet- 
ed as shown in Fig.1. Most of the power 
of the incident wave goes to the right hand 
waveguide, a small part is reflected back in 


Incidant TE 6; —s1_, Transmitted TE 
ain 


: (Vo. lo) 
Reflected TE, 
(YV Cora! for) 


y Scattered S 
y Backward TE “> cattered 
(Y Trea)! fnea)) Forward TE 
f (Voom), | (am)) 
cattered 
0 = Scattered 
Backward TM e@amaneean 


(V' am ie nm 
(nm),! (nm)) OW er syalleerayp! 
Fig. 1—Junction between a waveguide having a 


perfect flange and a waveguide having 
an oblique flange. 
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the same mode as the incident wave, and the 
remaining part is converted into many un- 
wanted modes, which propagate both forward 
and backward. In Fig.1, the letters in 
brackets are the mode voltages and currents, 
which will be defined later. 

The electric and magnetic fields are ex- 
panded into a series by the mode functions, ‘” 
namely: 


B= DV wee + di Ve) ext, 


HA,= Da Lay (2) Riciy ate ui Ta (z) hiti, 


C1) 
[Di di [iy BZ) Czciyy 


H,= x Via het, 


where E, and H; are the transverse com- 
ponents of the electric and the magnetic fieilds 
respectively; E, and H, are the longitudinal 
components of the electric and the magnetic 
fields respectively; Vi and J; are the voltage 
and the current of the mode respectively, 
and are designated as V;,; and J;,; for the 
case of the TE mode, and as V, and J) for 
the case of the TM mode. For the present, 
circular cylindrical coordinates (7, ¢, z) are 
used. e,; and h,; are mode vectors which 
satisfy the normalized orthogonality conditions 
and are given in the following expressions 
for the TM and TE modes respectively, 


C1) =i) X1z= Vz Ta, ”), 
(2) 
e=htyXt.= VT WM, ~) xi. 


where i. is the unit vector in the z direction. 
€ziy and hc) in equation (1) are given as 


9 
te eR LEC 
f= 1Wéq Pw, 


: (3) 
i tele 
[2] cae Thal, 


where the scalar functions Ta and Ty, are 
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ee =—VaJn x i sin 
@=PwaIn%X eae no, 
(4) 
Tra=atadn tan o> 
t= ta n¢ LI) cin MPs 


which satisfy the equation V,27;+%;7:=0, and 
Peay and qr; are normalizing factors. 

The components of the electric field of a 
mode are writen with Cartesian coordinates 
as (@,’, éy’,é.’), and with circular cylindrical 
coordinates as (e,’,é,’,e’). Then, the trans- 
verse electric field E,,’/ on the junction plane 
caused by the wave of a mode (V’,/’) is 


E’y={-—(V'e’, cos 6 + Ie’, sin 6)sin 
+V’e’, cos p}e*™ i, + Ai,, (5) 


in which the primes are put on the quanti- 
ties in the No. 1 waveguide (see Fig. 1), h is 
the propagation constant divided by the com- 
plex number 7, and i, and i, are the unit 
vectors in the ¢ and , directions respectively. 
In equation (5), the mode number 7 is omit- 
ted. If desired, the quantity A can easily 
be obtained, but since it is of no value in 
succeeding calculations, it is not necessary to 
find A explicitly. For a wave propagating 
in the +z direction, the upper part of the 
double sign in equation (5) is used and for 
a wave propagating in the —z direction, the 
lower part is used. 

The functions of 6 in equation (5) are ex- 
panded in power series of 0, and the terms 
of order 6°, 6! and 6° are collected giving: 


Blu=| V ey sin 20+ VeM +-dihr cos 0) 
oO” Oe 2 42 2 5 
~~ (sin g+h?rcos #5 


—Il'e',(+6—6 thr cos ¢)sin ¢ jis +Ai,. 


(6) 


In the similar way, the transverse magnetic 
field on the junction plane caused by the 


wave of the mode becomes 


Hy; Xig= 


-| rit | £1—dihreose# " + her eoste! 


+V'h',(6 +6" thr cos ~)cos ¢ jit Bi, (7) 


where B also need not be found explicitly. 

The electric and magnetic field components 
of each mode which is necessary in the suc- 
ceeding calculations are derived using e- 
quations (6) and (7). At first, for the TE; 
mode 


EB’ 3,9= Vo'arorteon{ Ji (X07) 


—dihtoyjr cos oJ, (%to11") 


52 
Zs 


{(sin’ o+h?[o137 cos’ ©) Ji (Xto117) 
+2 (or cos? gfi’(At0117)} jee 

11.) Xiz=— Vo'arorr%con{ — Yoo Ji A011) 
+6 | Yeontheonrh eon) 


Ke ; 52 

— 1 7 Xeon) |cos o+——Yrottour 
1 [lo 2 

xi! Xo) + Yond + 2? pour i Ato) 


x : - 
ee Cont 7, e137) jos? a Jet Ci, 


OLo 


(8) 


where C also need not be found explicitly. 
The normalizing factor g in equation (4) is 


Aue ik Pe 
nm]— ase Lan ? 9 
eiez / T V tnmy— 1? In Reman) 


where tym) is the mth root of Jn’(x)=0 for 


5o2 


all values of m except zero, Rinm]=*enmid, and 
a is the inner radius of the waveguide. The 
admittance Yt;) for the TE mode in equation 


(8) is given as 


pane A) (10) 
Hie Vii] Opto 
Secondly, the electric and magnetic fields 
on the junction plane caused by the reflected 
TEo; wave are 


E' a cou= V' c0199 019% 013 J: 20011) bo, 


HA’ ac) Xiz= — V’ cor 9t019% 0019 Y cor 1 A001) bo 
(11) 


In the same way, for the scattered backward 
TE, waves 


Enna — V'tumtind Ci Ji’ XtimY) (cos 


es ‘ Um] 
+6thtimyr cos? ©) +0Ytim}— 
2109 


xJiXmr) sin 2¢e Li, + Dir, 


Ja Opened Vcumadtind | ial 
xJi! XtimIY) (cos ¢ + dihtimr Cos? ”) 


X2 
+0 at) cos? O ! i,+ Ki,, 
1 [No s 
2) 


where D and E need not be found explicitly. 
In equation (11) and (12), it is taken into 
consideration that V’(o; is at most, the square 
of 6, and V’tm) is first order of 6.“ More- 
over, the amplitudes of the scattered back- 
ward TM;,; wave and the scattered waves 
other than TE,,, are zero,‘'? therefore these 
waves need not be taken into consideration. 

The continuity conditions of the electric 
and magnetic field at the junction plane are 
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E’ no +E’ a,co.+ 2 1,01mI1+ LE ni 
mn 
= Voertor + Ls Vier, 
4G 
CH 1149+ H n,c01 + DH a,timit+ 2D H'a,1) Xtz 
m } 


= ohio t+ DLjhis) Xt, 
J 
(13) 


where i represents any mode other than [01] 
and [lm], and 7 represents any mode other 
than [01]. The scalar product of equations 


(13) and groy%roJi(4ron) i, are integrated on 
the junction plane, giving 


V/A — Fo?) + WV’ fo + Wd 
m 
2 2ihrmyakcouk tim) _5— y, 
V Ft1mj—1 (Rpm —F*t011)? 


Vo’ Yt0111 — Fo?) — V’ p01 Y tor — 2 V tam) 
m 


7. OY tim] 


A tk om) 
heimya V Rt1mj—1 (Rm) — F?t01) 


x Yim} = Vo¥ ton, 
(14) 


where the F’s in both equations are equal. 
Therefore, they vanish in succeeding calcu- 
lations, and need not be calculated. Equations 
(14) are simultaneous algebraical equations, in 
which V’(o1] and V, are the unknown quanti- 
ties. By solving these equations, V’:o1] can 
easily be obtained. The reflection coefficient 
R, is the ratio of V’to1) to the voltage of the 
incident TE); wave. Using the mode con- 
version factors V’tim) from the incident TE, 
wave to the TE),, waves,‘? namely 


V 2 tkcoryk pm p01] — Atm) 2a } 
De. 2 2 2 OV’, 
2 VR tmi—1 Ram — R019) Atm) 


(15) 


Wm = 
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then the reflection coefficient becomes 


R,=Vto _ __ Prone 
We 2hio1a” 
Rum) 
xr ful — 
m (Ptimi—1) (Fm) —F?p019) Atm) 


(16) 


The reflection coefficient due to a slight tilt 
has been formulated as the summation of the 
functions of the eigenvalues of the modes 
generated at this slight tilt. In the calcu- 
lations of equations (5), (6) and (7), the 
following relations concerning the coordinate 
transformation have been used. 


C€=r(Cos ¢)9, 


Ppa r COS? g, 
2 G7) 
es 

g’=et+ 4 sin 29. 


The cross section of the No.1 waveguide 
with the junction plane is an ellipse and _ its 
area is a slightly larger than that of the 
No. 2 waveguide. But the area of the part 
of the No. 1 waveguide, not contained in the 
cross section of the No. 2 waveguide is only 
the order of 6”, and moreover the value gro1 
XrorJi1(“toy”) which is used to multiply the 
equation (13) to obtain the equation (14), is 
the order of 6 at least in the vicinity of the 
waveguide wall. Therefore, this area does not 
contribute of the order of 6” to the integral. 

We now show a numerical calculation for 
the reflection due to a slight tilt at the junction 
of the circular electric wave transmission line. 
For the case where the waveguide diameter 
is 51mm and the frequency is: 50 Gc/s, the 
terms following the sign >} in equation (16) 


m 
are divided into three parts and they are 
calculated separately as follows: 


Dee DL, areca (18) 


The first term concerns the modes which can 


propagate, and is a real number. The second 
and third terms concern modes which are 
below the cutoff frequency, and are imaginary 
numbers. The value m which is the bounda- 
ry between the second and the third terms 
does not have any special mathemetical mean- 
ing, but the larger m we choose, the higher 
the accuracy of the calculation becomes. We 
calculate the third term roughly as follows: 


ee sae: 1 
pe. ee 


a 


a ; 
as wii = E1505 X10", (19) 


in which the relation 


2 1 
Dg =} 20206, (20) 


m=1 


is used. After the calculation, the reflection 
coefficient is 


R.=(—0. 0259+ 0. 000268) 6”. (21) 


When the angle of the tilt is 10’, we obtain, 
using equation (21), a reflection coefficient of 
—133 dB. 


2. The Reflection due to a Slight Offset 


Reflection and mode conversion occur at 
the junction of a circular electric wave trans- 
mission line if the axes of the two wave- 
guides do not coincide. For this mode con- 
version, the first order solution had been 
obtained, moreover, it had been shown that 
the first order reflection is zero.“ Therefore, 
the reflection is much smaller than the mode 
conversion factor if the offset value is not 
large. But, in case of the helix waveguide, 
the reflection coefficient is an important value 
when we intend to investigate the delay 
characteristics of the waveguide, as explain- 
ed above. 

The electric and magnetic fields in the 
circular waveguide can be derived from two 
scalar functions Tium} and Tynmy. We will 
now transform the coordinates of the function 
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T in the No.1 waveguide as shown in Fig. 
2 into the coordinates of the No.2 waveguide. 


No.1 No.2 
Waveguide Waveguide 


Cie” ek ey 


Incident 


Fig. 2—Slight offset at a waveguide junction. 


Namely, from an additional theorem concern- 
ing the circular cylindrical wave,‘ 


: COs 
Thnm\=Qnm) Jn Aum) _ <2" 


= cos : 
= nm] ve Ji AtnmiD Jn+i Atami") ane (N+1)9, 
(22) 


where the primes are put on the quantities 
in the No. 1 waveguide. 

The fields in both waveguides are consider- 
ed as the superposition of all modes as shown 
in the lower right-hand corner of Fig. 1. 
Using equation (22) and the series expansion 


of the circular cylindrical function concerning 
d=a/a, 


Jihad) =Ji (Rid) =( “e yi ; 


ere! te ix seen ! 
(i+1)! ) ie 


=), il, 2 Mickeiercte ‘ 23) 


the fields of each mode in the incident wave- 
guide are expressed by the coordinates of the 
No. 2 waveguide, and are calculated up to 
the order of 6°. Using the definitions 


Fem? =3cam, 
(24) 
Fumio = 
2 =%nm], 
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we calculate the transverse electric field BE’. 
At first, for the incident TE); mode, 


E’' 11,0= Vo' aco (%c01 1 — Pt 1019) 
— 26po1%to11/1’ Aco) Cos Y 
+8 o13%t019 Jo’ (xcory7) cos 2p}i, } 29) 
+Gi,, 

Fl 11,0 X t= Y (011 E' 11,0, 


where G can easily be calculated, but need 
not be, because it is of no use in the suc- 
ceeding calculations. Secondly, for the re- 
flected TE); mode, 


E’ n,c01 = V’cor14corco1 Ji (A001) be, . 
(26) 
A n,001) X i2= — Yo 2", 0013. 


Thirdly, for the scattered backward TE)» 


modes, 
E’ am) = — V' tim giim imi thi’ Atm) Cos @ 
+6t1mJ1 Xt1mID +6 t1miJ2’ Xam”) 
x cos 2¢}i, + Hi,, | 


Ay ,01m) Xtz= — Yom’ atim), 
(27) 


where H need not be calculated also. The 
continuity conditions of the transverse electric 
and magnetic fields are 


FE’no+E’n,out+ DE’ n,m)+ DE a3 


m z 
= Voerory+ 2 V3e2;, 
J 
Wn 0+ Hn 0+ DH n,mit+ DH n,) Xiz 

m a 

= ohipout Shi) Xiz, 
Jj 
(28) 


where 7 represents any mode, except [01] 
and [1m], transmitted backward, and j re- 
presents any forward wave except [01]. 
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Both sides of the two equations of (28) are 
multiplied by — gporxcoui(xt0117) scalarly and 
integrated in the junction plane giving 


V,’ ar V’ ton] ie DRAG V' tm) = Wop 
m 


Yt011V0’ — YrouV’ to (29) 


te a! Cig lm V Tin) tO V oy 
m 
where 


Tin = 


V2 kro p19 
30 
V Rt1mj—1 (Roimi—R por) me 


Then the mode conversion factor from the 


incident TE), mode to~the TE;,, modes, 
namely 


Wty =(5 F [01] il ) 
a Yom) 


rok? 1mI0 Vi, 
a. 2 VRom—1) pm) — p01) 


G1) 


which is given in reference,‘ is substituted 
into equation (29), and equation (29) becomes 
two simultaneous equations having two un- 
knowns, V’(o1j and V)’. We can easily obtain 
V'tory by solving these two equations, and 
the reflection coefficient R, is the ratio of 
V'to1y to Vo’, namely 


Re V'co1j_ tor”. tor re) 
Vi’ o m heim] Aron) 


eras 
x 3 RY, 
(R°t1m]—1) (tm — R201)? of 


The reflection coefficient has been obtained 
in the summation form of the functional 
terms of all eigenvalues of the modes which 
are generated by this offset in the first order. 

We now proceed to calculate numerically 
using equation (32) for the case where the 
waveguide diameter is 51mm and the fre- 
quency is 50Gc/s. The term following the 
sign >| in the equation (32) is divided into 
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three parts, and they are calculated separate- 
ly as follows: 
8 S 
DS) see Te (33) 
m m= m=17 
The first term concerns the modes which can 
propagate and is a real number, while the 
second and the third terms concern the 
modes which are below the cutoff frequency 
and are therefore imaginary numbers. The 
value m which is the boundary between the 
second and the third terms does not have 
any special mathematical meaning, but the 
larger m we choose, the higher the accuracy 
of the calculation becomes. ‘The third term 
does not converge if the following rough 


estimation is used. That is, for large m, 
using the relation krj,_j)mz, we have 


ay ee | 1 
BETO an a 


m=17 m= hoa Rom) 


= Dee) — Divergent. (34) 
heojaz m=17 M 


The divergence occurs because only the first 
term of the power series (23) is used at the 
time of the coordinate transformation, how- 
ever large k; is. To prevent this divergence, 
the third term is calculated from the original 
cylindrical functions as follows: 


ss B 4 Z Ji:(mn0) 2 
2-15 yy 
may hpovjax m=17 =m eS 


Since 6X1, the right-hand side of equation 
(35) can be obtained using the following 


integral relation: 


Sy Lim ag” Lae 
pro 


n=pe —H0 ne 
Sd aye log (px0)Jo(pr0) 
2pro 2 
lye ‘ 
— 5 |" los xtra) (36) 
Z pro 


If pzd<1 is assumed, the integration of the 
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last term of equation (36) can be performed 
changing the lower limit, p70, to zero, with- 
out serious error, because the integrand is 
very small in the vicinity of s=0. The inte- 
gration in such a case is 


|; 108 i iCNdee esos 2 (37) 
0 


where 7 is Euler’s constant, namely 7= 
0.557721. Using the integration relation, 
equation (36) becomes 


gy Film) «9 604—0. 52 log (pxd). (38) 
mo 


When the frequency is 50 Gc/s and the wave- 
guide diameter is 51 mm, 


Ro= (0. 126—7{0. 154—0. 0888 log (1770)} Jo”. 
(39) 


If the offset is 0.1 mm, for example, the re- 
flection coefficient is —108 dB. 

When the mode conversion factor (31) is 
calculated,“” only the first term of equation 
(23) is used. Therefore, for very high order 
modes, large errors occur: 


3. Mode Conversion at the Junction be- 
tween Circular and_ Elliptically 
Deformed Waveguides 


Mode conversion occurs at the junction 
between a circular waveguide and an eellipti- 
cally deformed waveguide having the same 
peripheral length as the circular waveguide, 
as shown in Fig.3. There are many types 


; No.1 No.2 
Incident Waveguide Waveguide 


Fig. 3—Elliptic deformation discontinuity 
at a waveguide. 
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of deformation of the cross section of the 
waveguide, but elliptical deformation is most 
liable to happen in practice and is easiest to 


analyze. 
The 7 functions of each mode in the 


oe : : ; 
elliptic waveguide are given as‘” 


T’t01] =4' on Jo’ tour’) 


12, 


(40) 
5 ~ COS 
T' (nm) = Ol tnmi In tami” ) —sin ng’, 


= sin 
T’ nmy =D! cnmy In ® nm”) Fae no’, 


where the primes are put on the quantities 
of the No. 1 waveguide, and e is the elliptici- 
ty. In equation (40), the necessary terms to 
calculate the mode conversion factor up to 
the order of 6 are used. The J’ functions in 
the No. 2 waveguide are 


Tro = 40011 Jo“ 011”), 


4 cos 
T (nm1=Qinm) Jn Ctam”) ere 
= Silis (41) 


, sin 
des = News Cera No. 
cos 


From equation (2), the continuity conditions 
of the electric and magnetic fields at the 
junction plane become 


Vo'ViT’ (01) Xiz+ VrorVil’ X22 
+21 temiVel famixts oo ee eee 
=VoVeT 01] Xie + DO VinmVeT tam) X tz 
+ i VamyViT cam, 

Yt Vo Ve T” ora X be — Y’ cor tore T” t01 X iz 
LV tamiV tam iT’ tam) 


=F Da YF dias YW Guns Ve Le (nm) = Yto1) VoV: T{01] 


a ee 
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Xtz+ Di Ytum) VinmVt T lnm] X tz | 
ai Ds Vaams Vian y IIPS. 
(42) 


The electric and magnetic fields in both wave- 
guides are represented as the superposition 
of all modes as shown in the lower right-hand 
corner of Fig.1. The value A in Fig.3 is 
related to the ellipticity in the equation 4A/a 
=e", and we define 6 by 6=2A/a. At first, 
the first equation of (42) is multiplied by 
OtamV i {J2(Xt2mi7)cos2¢} Xi, scalarly and is inte- 
grated in the junction plane. Neglecting the 
terms of order higher than 6°, we get 


Rk’? . P 
V0'q' torr ee acona\ \¥i {Jo(X' tory cos2g} x iz 


Vi {Je (X(2m]7) COS 20} Xi.dS+ V’ tam] = Vom. 
(43) 


If the same operations are performed on the 
lower equation as on the upper one, we 
obtain equation (43) in which the first term 
of the left-hand side is multiplied by Y’fo1, 
the second term is multiplied by Y’~2m3, and 
the right-hand side ‘is multiplied by Yf2mj. 
By integration in the region which is not 
common to the circular and elliptical cross 
sections, we do not get the term 0(0), there- 
fore the integration range of (43) may be set 
equal to the circular cross section. From 
equation (43) and the result calculated from 
the lower equation of (42), we get 


Ve2mj— V' t2m]= Ln Vo' e’, 
(44) 


Voom tom t Ym tami= Le’ to V0 &*, 


where 


ee V2 Rrosk tem) (45) 
A / Fam] —4 (Fam — R01) 


Equations (44) are simultaneous algebraic 
equations having two unknowns, Vj2_) and 
V’t2m}, and the solutions of these equations 
are 


Y’ tam] + Y’ ton 
View = jie Len V5 C2 
Bol Vine) tol ge 
(46) 
Yu) (ol 4 
Vom] = dimV oe. 
ISN STEN Ge ae 


The mode conversion factor from the incident 
TE); mode to the scattered forward TE», 
mode, Ct); and the mode conversion factor 
to the scattered backward TE2,, mode, C-t2m7; 
are given in one equation 


Agl2mj= Aglo1) 
2 v7 Aglorol2m) 


(Cama = Mop @ (47) 


where Ln» is a coefficient not dependent on 
frequency nor dimensions and is shown in 
Table 1. In the case where the frequency is 
50 Gc/s and the waveguide diameter is 51 mm, 
and when there are discontinuities of 2A= 
0.1mm at intervals of 5m, and if the con- 
verted power is immediately transformed into 
heat loss due to the conductivity of the wall, 
then the mode conversion loss is 0.065 dB/km. 
As is easily calculated from equation (40), the 
mode conversion to other than the TE»,, 
mode does not occur in the first order of 0. 
And also, reflection does not occur in the 
first order of 0. 


Table 1 


Lm IN EQUATION (47) 


Mode | TEs; 


TExs | TEs, | TE, | TEx 


| TEx | TEx: i, 
| | 


bee —1,024 0. 313 a 0. 163 


— ; i 
0.114 | 0. 088 | 0. 071 | 0. 062 


| 0. 054 
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4. Mode Conversion in the Vicinity of 
the Cutoff Frequency 


As the frequency approaches the cutoff fre- 
quency of some mode, the mode conversion 
factor increases infinitely for any case of 
mode conversion. This corresponds to the 
phenomenon in a cavity where the field 
strength increases infinitely as the frequency 
approaches the resonant frequency if the 
dissipation loss is neglected. It may be con- 
sidered that a waveguide at the cutoff fre- 
quency is a resonant cavity in which the 
resonance is independent of cavity length. 
We will now calculate the exact mode con- 
version, taking the dissipation loss into con- 
sideration, when the frequency is near the 
cutoff of some mode, for the case of an 
elliptically deformed discontinuity. 

The mode conversion factor has been shown 
in equation (47), with the assumption of zero 
dissipation, but when dissipation is taken into 
consideration, this factor becomes 


; Lime? hoot hem) 
2 / hoi} V him) 


C*t2m]= (48) 


From equation (47), the mode conversion 
factor when the frequency is much _ higher 
than the cutoff frequency is 


Com iar = 1 me", (49) 


and from equation (48), when the frequency 
is near the cutoff frequency, it becomes 


Cin Ce a (50) 

V htom] 
The eigenvalues of Ly of the equation (45) 
become complex numbers if the dissipation 
loss is taken into consideration, but ZL,» for 
the case of no dissipation has been used be- 
cause the imaginary parts are negligible 
compared with the real parts. The ratio of 
the mode conversion factor near the cutoff 
frequency to the value at frequencies much 
higher than the cutoff frequency is 
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'C* tami — 


=i Talia ee 
yen: C [2m] ai ee (a?(2m] + P[2m)) a 
[2m ]far 


C* (2m]rar 
1) 


where a2m) and ft2m] are the attenuation 
constant and the phase constant of the TE2,, 
mode respectively, and when the frequency 
is near the cutoff frequency, they become‘ 


Qinm] = y= @+S)+ Jee 


j (Gy) 
een (ES ya Meee . | 
where x and S are 
x= hp 2S fo ) 
(93) 


2 we f \? n2 } 
S=—— 1 ( 
6g a ea, aE Rtnmj—n* {? 


and where f, is the cutoff frequency; & is 
the phase constant in free space; 6 is the 
skin depth; ¢o is the conductivity; and a is 
the waveguide diameter. Then T in equation 
(51) becomes 


Pee een {GE Se Se (54) 
This value becomes maximum when x= —S, 
and it is 

Ue at Sissies (Gey 


When the waveguide diameter is 51 mm, the 
cutoff frequency of the TE, mode is slightly 
lower than 50Gc/s. At a frequency a little 
lower than the cutoff of the TEs mode, 
namely when the relation «=—S is satisfied, 
T becomes maximum and is 8.5 for the case 
when the relative conductivity is 100 %. The 
mode conversion loss was 0.065 dB/km_ for 
the elliptically deformed discontinuity of 0.1 
mm at intervals of 5m as shown above. At 
the frequency at which T becomes its maxi- 
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mum value, the mode conversion to the 
forward and backward TEss waves increases, 
therefore the mode conversion loss increases 
by 36 % and becomes 0.088 dB/km. 

Here we studied the effects of the cutoff 
frequency only for the elliptically deformed 
discontinuity, but similar phenomena would 
occur for the case of other types of discon- 
tinuities, and similar analyses could be made. 


Conclusion 


A method to calculate the reflection and 
the mode conversion due to a slight discon- 
tinuity at the junction of circular waveguide 
has been derived, and the reflection due to 
a slight tilt and a slight offset, and the mode 
conversion due to a slight elliptically deform- 
ed discontinuity have been obtained. It has 
been proved that the reflection coefficients at 
the slight discontinuity are very small in the 
above two cases, but these must not be 
ignored when the delay distortion of the 
circular waveguide transmission line is analyz- 
ed, especially in the case of a helix wave- 
guide. fy 

The calculation of the mode conversion due 
to a slight elliptically deformed discontinuity 
would be useful to determine the necessary 
accuracy of the dimensions and_ installation 
of the waveguide transmission line. The 
mode conversion loss was 0.065 dB/km, due 


to the discontinuity of 0.1mm _ successively 
located at an interval of 5m, and at the 
frequency near the cutoff of the TEs, mode, 
it increased to 0.088 dB/km. 
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Matching between Traveling-wave Tube Helix 


and Waveguide’ 


Tomoyuki UNOTOROt 


The method of orthogonal array layout of experiments was applied to the determination 
of dimensions of the three factors which can be ‘chosen independently of the required 
performance of a traveling wave tube and affect the matching between the waveguide 
and the helix. As a result, these factors were found to affect the matching only inde- 


pendently. 


A broad band matching has resulted from the design thus made over the 


frequency range of up to 20 percent of 11 Gc/s. In addition, the contribution of each 
wide factor to the admittance cemponents has definitely been determined. 


1. Introduction 


The helix type traveling-wave tube has 
broadband characteristics. The practical fre- 
quency band, however, is limited by the 
matching between the traveling-wave tube 
helix circuit and the input and output coupling 
circuits. The matching method may be divided 
into two main categories: one in which a 
coupling helix is used for connection to a 
coaxial circuit, and the other in which a short 
antenna is provided at the end of the helix 
and inserted in a waveguide. The matching 
method using a coupling helix, which results 
in an extremely broad bandwidth, i. e. a band- 
width of an octave or more, is widely used 
in measuring instruments and for other 
purposes. However, since some residual _re- 
flection due to mismatching occurs in the 
pass band, this matching method may not be 
used for multichannel communications systems 
with an extremely large number of channels, 
in which even slight phase distortion will 
cause cross-talk troubles. The second method, 


* MS in Japanese received by the Electrical Communi- 
cation Laboratory, April 18, 1960. Originally publish- 
ed in the Kenkya Zituydka Hokoku (Electrical Com- 
munication Laboratory Technical Journal), N.T.T., 
Vol. 9. No. 6, pp. 639-648, 1960. 

| Electronic Components Research Section. 


hh 


whereby connection is made to a waveguide, 
will provide a somewhat narrower bandwidth, 
but as it is possible to make the residual re- 
flection in the pass band very small, this 
method has been adopted for the most part 
in the ultra multichannel microwave relay 
systems now in general use. The coupling 
structure connecting the helix and the wave- 
guide has thus far been empirically designed. 
Therefore, it has not necessarily been possible 
to determine whether a design was optimum. 

In this study, the factors controlling the 
matching characteristics, have been qualitative- 
ly ascertained as a result of systematic ex- 
periments on the coupling circuit between the 
helix and the traveling-wave tube as a_ part 
of 11 Gc band traveling-wave tube develop- 
ment researches. These results are applicable 
not only to 11 Ge tubes but also to tubes for 
other frequencies, as long as they have the 
same structure. The conclusions should also 
prove valuable in the design of coupling 
circuits. 

In this coupling circuit, as described in 
detail later, a short antenna is connected to 
the end of the helix circuit and the other end 
of this antenna is connected to a_ small 
cylinder. The circuit is inserted into the 
center of the waveguide. That the conduct- 
ance component of the coupling point ad- 


VOLUME 8, NUMBERS 11-12, NOVEMBER-DECEMBER, 1960 561 


mittance varies with antenna length, while 
the susceptance component of the coupling 
point admittance varies with cylinder length. 
Each component is independent, there is no 
interaction between them. 

A coupling circuit was thus made with the 
best dimensions obtained by making use of 
this result, so that good  ultra-broadband 
matching characteristics could be obtained in 
the 11 Gc band over a bandwidth of 2200 
Mc. It could also be deduced by the investi- 
gation made based on the method of. statis- 
tical analysis of an orthogonal array layout 
of experiments that matching giving a volt- 
age standing wave ratio of less than 1.3 could 
be obtained without adjustment in mass pro- 
duced tubes with less than 2 % loss. 


2. Coupling Structure 


Many methods have been used for coupling 
between a traveling wave tube helix and a 
waveguide. One of the structures much used 
is a tapered helix,“ in which the end of 
the helix is straightened gradually into a 
straight wire which is inserted into the wave 
guide. Another is a straight antenna,‘ in 
which a short antenna is provided directly 
at the end of the helix and inserted into the 
waveguide. In the former type, it is difficult 
to determine the dimensions of the tapered 
part where end of the helix is to be straight- 
ened gradually, and furthermore, this part 
does not contribute to the gain of the tube, 


Waveguide 


Glass rod 


 —— at 


CMG 
Shield pipe 


Short 
plunger 


Fig. 1—Coupling structure between 
helix and waveguide. 


but instead has a harmful effect on the focus- 
ing of the electron beam. And, moreover, 
when tape is used for the helix, it is difficult 
to straighten. For this reason, the latter type 
was used in this Laboratory for all the tubes 
that were developed. 

The structure of the straight antenna type 
coupling has been designed, taking has been 
into account mechanical requirements, as 
shown in Fig. 1. 

The helix circuit is supported by three 
quartz or ceramic rods, whose ends are insert- 
ed into three holes made in the choke. The 
helix end is connected directly to the antenna 
running parallel to the tube axis, and the 
antenna is connected to a hollow cylinder 
jutting out from the choke. This assembly 
is inserted in the center of the waveguide, 
as shown in the figure. One end of the wave- 
guide is shorted with a short plunger, so that 
microwave power comes in or out through 
the other end. 

The radius of helix is determined by the 
electronic design of the traveling-wave tube, 
the diameter of glass tube is determined by 
mechanical dimensions of supporting insulator 
rods, and the height of waveguide is limited 
by the focusing magnet. Therefore, there 
are only three coupling circuit dimensions 
which can be freely changed, without regard 
to these requirements, for the purpose of 
controlling the microwave matching charac- 
teristics. They are antenna length, cylinder 
length and position of the short waveguide 
plunger. In some cases, a tuning mechanism 
is provided at waveguide circuit outside the 
tube, thus making it possible to effect match- 
ing by adjusting the tuning. This method, 
however, is not desirable, since it will es- 
sentially deteriorate frequency characteristics. 
The consideration, therefore, should be to 
make out the design of the coupling circuit 
without reflection. Systematic study has been 
thus made as regards above mentioned three 
dimensions, which can be changed in the 
tube by itself, on the standpoint such as how 
they affect matching characteristics and, further 
whether complete matching can be obtained 
by merely adjusting these dimensions. 

The samples measured are intended for 11 
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Table 1 


DESIGN OF 11 Gc TW... CIRCUIT (mm) 


Tape helix dimensions (hae Wh 
Mean helix radius 253 
Helix pitch 0.9 
Inner diameter of shield pipe 10. 0 
Coupling waveguide dimension 5.2X22.9 


Gc band high power tubes. The dimensions 
determined in designing traveling-wave tube 
are as shown in Tablel; while cylinder 
length J, antenna length g, and the distance 
from the glass wall of the tube to the short 
plunger p shown in Fig. 2 represent the di- 
mensions that can be changed. Since the 
diameter of the glass tube is 8.5mm, adding 
4.25mm to the value of p will give the 
distance from the center of the tube to the 
short plunger. 


beer 
| 
_ fH 


Fig. 2—Coupler dimensions 
to be changed. 


In making this study, not only admittance 
was sought for, but also statistical analysis 
of the orthogonal array layout of experiments 
was made. The question, therefore, was how 
many levels should be chosen for the above- 
mentioned three dimensions. If only two 
levels were taken, only the first order tenden-: 
cy of effects could be obtained. In this case 
it would not be possible to determine the 
optimum level, although it would he possible 
to determine a good one. In order to obtain 
the optimum value, at least three levels are 
required, so that the second order appro- 
ximation can be obtained. The increase in 
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number of levels will raise the accuracy. 
However, if four levels are used, 64 sets of 
experiments will be required, so that a great 
deal of labor will be necessary. 

It was, therefore, decided to make experi- 
ments with three levels. In choosing these 
three levels, it would be meaningless to take 
them far from the point which gives good 
matching. In our case, however, it was not 
known at first which point would give good 
matching, so that levels over a rather wide 
range were used. The values to be adopted 
for the analysis of variance of the orthogonal 
array layout of experiments were determined 
after the results of the measurements were 
examined. 

The design values of samples and the values 
obtained by measuring the dimensions of 
samples are given in Table 2. In manufactur- 
ing the samples, an ordinary industrial pro- 
duction process was followed, without any 
special attention and the order of manu- 
facturing was random. 


Table 2 


DESIGN AND MEASURED DIMENSIONS 
OF SAMPLES (mm) 


anor | Design Measured 

l g l g 
0.5 ate 0. 66 2.44 
B 0.5 <5 0) 0. 54 3. 56 
C | 0.5 4.5 0. 62 4. 62 
D We 15 | 1.54 92. 
E eS PAIS Mt 2. 89 
EF is 5) 35S) 1.49 3. 36 
G Bi) 4.5 1554 4. 68 
H | eo) OFS ZaOe, 0. 67 
[i 2, Meh aanee 2.45 1.68 
AI 2.5 Fah 2,03 2a ie 
K Be) B55) 7o5 (33 3. 58 
lie, BES 4.5 2.61 4.59 
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The position of inserting this coupling into 
waveguide should be so arranged that, as 
seen from the waveguide, the upper face of 
the choke at the joint of the cylinder is on 
the same horizontal plane as the inner wall 
of the waveguide. Although errors are un- 
avoidable in doing this, the statistical study 
is industrially very effective, as the random- 
ness of the results due to such errors can 


thereby be deduced. 
3. Measurements 


In making the measurements, the position 
of the short plunger in each sample changed 
by 1mm at a time, and the admittance 
measured from the waveguide side, with the 
frequency changed markedly in each case. 
Since the helix circuit is very heavily loaded 
it may be considered that a virtually reflec- 
tionless termination is connected to it. 

The measurement assembly consists of a 
coupling waveguide with a height of 5.2mm, 
which was, as shown in Fig. 3, connected 


Reference short 


Twisted waveguide 


Direct 
admittance 


Waveguide switch 


measurement 
set length 40. 5mm 


Taperedeuce= al 


{Short plunger 


Fig. 3—Measurement assembly. 


through a tapered waveguide 40.5 mm _ long, 
to the 11 Ge band standard waveguide, i. e. 
10.2 X 22.9mm waveguide. No tuning mecha- 
nism was used so as not to cause any reflec- 
tions. The measured admittance was _trans- 
formed into the value at the coupling point 
on the waveguide, thereby making it possible 
to examine the effect of the change in di- 
mensions. In making statistical analysis of 
the experiments, only the voltage standing 
wave ratio was examined. 


4. Analysis of Variance and Estimated 
Value by Orthogonal Array Layout 
of Experiments 


The method of orthogonal array layout of 
experiments is a very effective means by 
which to obtain the combination of factors 
to give the optimum results by averaging the 
random errors of all factors, though there is 
anavoidable randomness in many factors that 
are considered to affect the final product in 
the production process. 

In the case of this traveling wave tube 
coupling circuit, it is important not only to 
experimentally manufacture it carefully in the 
laboratory with an eye to obtaining good 
characteristics, but also to have it produced 
in the factory; where it will be necessary to 
cope with uncertain factors which might arise 
at any time and which are beyond the con- 
trol of a mere gauge or jig. As the statistical 
estimation is considered quite effective in such 
a case, the following examination was con- 
ducted by making use of this method. 

It is necessary, in making analysis, to ex- 
press the matching quality numerically for 
each sample. Although it is indisputable that 
the quality of matching in each specific case 
may be indicated by the voltage standing 
wave ratio, opinion may be divided as to 
how the quality of matching of a certain 
sample for the whole frequency band should 
be expressed numerically. Here, however, the 
following tentative plan is proposed as a scale 
for expressing such quality. Namely, the 
standing wave ratio is taken at about 18 points, 
with the frequency changed within the band 
from 10.6 to 11.44 Gc., and the values of the 
standing wave ratio at the various points are 
averaged. If, in this case, a simple mathe- 
matical mean is taken, it will indicate a good 
value as a result of taking an average as a 
whole, even if there are considerably bad 
values at one or two points. Therefore, so 
the value thus obtained will not coincide with 
the quality determined by the observation of 
the measured points on an admittance chart. 
However, by taking a geometric mean, a 
reasonable value could be obtained which 
coincided well with the determination by an 
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admittance chart. The geometric mean, there- 
fore, was adopted. To state the calculation 
method concretely, a mathematical mean 1s 
taken of the values obtained by squaring the 
values of the standing wave ratio (>1) minus 
i, and, with the square root of such mathe- 
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matical mean extracted, 1 is added to the 
value thus obtained to make up a weighted 
geometric mean of the standing wave ratio. 
Needless to say, the smaller the value thus 
obtained, the better the matching will be. 


Table 3 


SQUARED MEANS OF STANDING WAVE RATIOS 


Short plunger position (mm) 


Sample l g 0 1 2 a 4 5 

0.5 2.9 1. 662 1. 598 1. 918 DA SG 
0.5 3.9 1.814 1. 443 1.319 1,273 1. 818 

Cc 0.5 4.5 2. 055 1. 566 1. 602 

D 1.5 15 1. 455 ih, oy? 1. 398 

E is) 250 1. 334 1. 412 1. 207 1. 388 

FE 159 SK) 1. 349 1.159 1. 265 

G 11s '5} 4.5 1. 946 1.511 1.411 1.610 

ist eS) 0.5 2. 690 3. 450 4.118 

I 2.9 1.5 1. 432 LOLS 1. 750 

J 7) Ds 1. 353 1. 206 PeroT 

K Dh, 83 3.95 1. 238 1.135 1, 222 

L Bo 4.5 1.475 1. 261 1. 412 


Table 3 shows the numerical values obtain- 
ed by calculating the experimental results in 
this way, where the indicated position of the 
short plunger, as mentioned before, is the 
distance measured from the glass wall of the 
tube. As in the analysis of variance, only 3 
levels were taken; 3 positions of the short 
plunger sufficed. Here, however, many positions 
were taken in some cases, in order to find 
out the best point for matching. As it is not 
necessary for analyzing variance to fix the 
level of each factor, such level may be select- 
ed at some point in close proximity to the 
best point. The values to be adopted for 
analysis may be determined by examining 


the results shown in this Table. 

Table 4, which shows these numerical 
values layed out on the three level orthogonal 
array,*? comes next. Samples A, B,C...... 
have the dimensions given in Table 3, where 
the small suffixes indicate the positions of 
the short plunger ». As mentioned above, 
the values in Table 3 which are in close pro- 
ximity to the best point, namely the bold- 
faced numerical values in the table have been 
adopted, whereas in Table 4 such numerical 
values minus 1 are used for the sake of con- 
venience of calculations to be made later. 
Table 4 represents the 3 level 27 experiments 
type orthogonal array, where the experi- 
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Table 4 
LAYOUT OF EXPERIMENTS ON THE “L24(3!8)” TyPE ORTHOGONAL 
ARRAY AND MEASURED DATA 
No. of l g US eres Pp lxp gxp Measured Squared 
experiment 1 2 3 5 6 a Bo veeeee Sample mean value value 
1 1 1 1 1 2 1 Ay 0. 662 0. 4380 
2 1 1 x 2 2 2 AG 0. 598 0. 3570 
3 1 1 1 3 3 3 WW Se cAY 0. 918 0. 8440 
4 1 2 2 1 1 2 i es 0. 319 0. 1015 
5 1 a) 2 2 3 By 0. 273 0.0745 
6 1 2 2 3 3 1 B; 0. 818 0. 6680 
7 1 3 3 1 1 3 G3 1. 055 1. 1120 
8 1 3 3 2 2 1 ~Cy 0. 566 0. 3200 
9 1 3 3 3 3 2 CG 0. 602 0. 3620 
10 2 if 2 1 2 1 E, 0. 412 0. 1700 
11 2 1 2 2 3 2 je) 38 0. 207 0. 0428 
12 2 1 2 3 1 3 Ne aes 0. 388 0. 1506 
13 D ee 1 2 2 F, 0. 349 0. 1220 
14 2 2 3 2 3 3 F, 0. 159 0. 0253 
15 2 2 3 3 1 1 F; 0. 265 0. 0700 
16 2 3 1 1 2 3 mee. 0. 511 0. 2610 
17 2 3 i 2 3 1 Gi 0. 411 0. 1688 
18 2 3 1 2 2 Gs 0. 610 QO. 3787 
19 3 1 3 if 3 1 aad: 0. 353 0. 1249 
20 3 1 3 D 1 2 an F 0. 206 0. 0426 
ml 3 1 3 3 2 3 i 0.191 0. 0366 
22 3 2 1 1 3 2 Ke 0. 238 0. 0569 
23 3 2 1 2 1 3 ie 0.135 0. 0183 
24 3 2 1 3 2 1 K; 0. 222 0. 0487 
25 3 3 2 il 3 3 iis 0. 475 0. 2252 
26 3 3 2 2 1 1 ies 0. 261 0. 0682 
Di 3 3 2 3 2 2 ip, 0. 412 0. 1698 
| ‘Total Tis6l6 6. 4524 
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Fig. 4—Effects of each factor on the average standing wave ratio. 


Table 5 


ANALYSIS OF VARIANCE 


| Sum of 5 
Factor Freedom Puaces Variance F test 
1 Y 0. 66 0. 33 15, Ax 
g 2 0. 25 0. 125 beat 
p 2 0. 18 0. 09 4, 20* 
lxg 2 0. 03 0. 015 0. 70 
lxp 2 0. 01 0. 005 0. 23 
exp 2 0. 01 0. 005 0. 23 
Error 14 0. 30 0. 021 
Total | 26 1. 44 


1% level of significance. 


ments are arranged so that in columns 1, 2, 
and 5, the main effect of J, g, and p; and in 
columns 3, 6, and 8, the interaction of J and 
g, land p, and g and p respectively might 
be introduced. 

The analysis of variance was made _ using 
this table, and the finally obtained is shown 
in Table5, with the process of calculation 
omitted. The most important conclusion 
shown by this table may be derived from 
the results of F test, which mean that the 
effect of each independent factor is great, 
but no interaction is seen between any of 
the factors. In other words, it may be said 
that the effect of each factor is independent. 


It is, of deep significance that such a con- 
clusion has been obtained on a sstatistical 
basis. 

Fig 4 shows the effect of each factor graphi- 
cally. The reliability limit determined separate- 
ly, by making use of “Table t”, at 0.145 
with 99 % confidence and at 0.086 with 90 % 
confidence. It may also be seen from this 
graph that the difference in each dimension 
by 0.1mm or so from the best point will 
have only a small effect on the standing 
wave ratio. 

When sample K, which turned out to be 
the best as a result of the above-mentioned 
analysis, i.e. with /=3, e¢=2 andy f)— 7a 
adopted, the values with which the standing 
wave ratio is statistically expected to come 
out will be: 


Central value 1.039 
Upper limit with 99 % confidence 1.311 
Upper limit with 98 % confidence 1.279 


It was confirmed that from these numerical 
values, the weighted geometric mean of stand- 
ing wave ratios of below 1.3 can be obtained 
with 98% confidence, taking into account 
the randomness in the manufacturing. It may 
be seen also from Fig. 4 that the sample is 
quite promissing for industrial production, as 
the effect of errors in each dimension does 
not seem to be so significant. 

Since the above deduction is entirely based 
on the weighted geometric mean value, it 
should not be made a direct conclusion, but 
should rather be considered a statistical proof. 
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From this study, however, it has been made 
clear what combination of three dimensions 
will give the best matching statistically. By 
only such a method as described in the follow- 
ing section, whereby the admittance value is 
to be estimated, it would be difficult to de- 
termine the best combination mentioned 
above, since a distinct tendency of the effect 
of each factor as given in Fig. 4 may not be 
obtained. It may be said that this, together 
with the conclusion of non-interaction, is the 
significant result of the statistical analysis of 
the orthogonal array layout of experiments. 


5. Analysis by Admittance Chart 


In the preceding section has been described 
how the measured values of voltage standing 
wave ratio were dealt with statistically. 
Here, the effect of the dimensions of each 
part is investigated by entering in a Smith 
chart the admittance at the waveguide helix 
coupling. As for the admittance obtained by 
actual measurement, the positions of points 
on the chart are in disorder because of the 
effect of very small random reflection at the 
helix, in which case the better the matching, 
the greater such effect will be. Therefore, 
in finding the tendency of admittance variat- 
ions due to changes in dimensions, it would 
be impossible to grasp the general effect of 
such variations, only by looking these charts. 
In order to see the tendencies clearly, the 
admittance at 3 points; i.e., at both ends 
and the center of the frequency band to be 
measured were examined. 

First, admittance variations due to shift of 
the position of the short plunger were ex- 
amined. On Fig.5 are plotted the effect 
caused by the position of the short plunger 
as regards two samples A and H. Points 
with the numeral 1, 2, or 3 in the chart 
indicate the admittance at the coupling in 
the case of the frequencies of 11.4 Gc, 11.0 
Ge and 10.6 Gc, respectively. Small suffixes 
given to sample symbols indicate, as hereto- 
fore, the positions of the short plunger p. 
As, for instance, means Sample A with di- 
mensions of g=2.5 mm, /=0.5 mm, and short 
plunger position 2mm from the glass wall of 
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Fig. 5—Change in coupling point admittance 
depending on short plunger position. 


the tube. 

As seen from Fig. 5, it has been confirmed 
that the admittance at the coupling is not 
only loaded in parallel with the susceptance 
produced at the short plunger, but that the 
two factors of conductance and susceptance 
are affected by the position of short plunger. 
Since the diameter of glass tube is, as mention- 
ed before, 8.5mm, the distance from the 
coupling point, in case the position of the 
short plunger is about 4.5mm, will be 1/4 
wave length, so that the parallel admittance 
produced on the side of the short plunger 
will be infinitely great. It should, however, 
be noted that the point to give good match- 
ing is, as seen from Fig. 4 as well as from 
many admittance charts, nearer than 1/4 
wave length. 

Then, investigations were made as to how 
the admittance produced at the coupling point 
varies when other factors, i.e. antenna length 
g and cylinder length J are changed. Fig. 6 
shows, as in Fig. 5, a Smith chart on which 
the admittance measured at 3 points in fre- 
quency band for each sample is plotted. 
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Fig. 6—Change in coupling point admittance 
due to change in antenna or cylin- 
der length. 


Here, however, having an eye only to the 
‘effect of antenna and cylinder length, the 
position of the short plunger is fixed at 3mm. 

A very noteworthy characteristics is clari- 
fied from Fig.6. Namely, as seen from this 
figure, the variations of admittance at the 
coupling point when antenna length g is 
changed and when cylinder length Z is chang- 
ed are orthogonal to each other on the Smith 
chart. For the convenience of comparison, 
the layout of g and J for each sample is 
shown in Table 6. It may be seen by com- 
paring this table with Fig.6 that the arrange- 
ment of sample dimensions completely coin- 
cides with that of the admittance points. 
This fact is also supported by the conclusion 
reached previously through the method of 
statistical analysis of the orthogonal layout 
that there is no interaction between the 
effects of dimensions. 

Now that the change in antenna length 
mainly affects the conductance part, while 
the change in cylinder length mainly affects 
the susceptance part, the effects of the an- 
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Table 6 
LENGTH OF ANTENNA AND CYLINDER SAMPLES 


—~__ Cylinder length 


ret Dabs 1.5 0.5 
Antenna length —~|__ ———— 
4.5 eae ® G4) C 
Se —— 
3.5 K E B 
2.5 J 1s A 
ie 
Woes I D | 
0.5 | # | 


tenna and the cylinder at the coupling point 
may be understood as follows. Considering 
the coupling structure for the helix and the 
waveguide as an equivalent transformer, the 
change in antenna length will correspond to 
a change in the ratio of transformer, and 
thus the residual inductance due to leakage 
flux of the transformer will be neutralized 
through the capacitance introduced at the 
cylinder. In fact, it is well known that if a 
cylindrical post is provided inside the wave- 
guide, it will introduce capacitance. There- 
fore, although the above-mentioned experiment 
made with respect only to one dimension of 
helix and waveguide, it may well be conclud- 
ed that the obtained effect of- antenna and 
cylinder length on the admittance may 
generally exist, regardless of the dimensions 
of helix and waveguide. 

A similar conclusion has been reported” 
with respect to a tapered helix type coupling 
circuit, in which case, however, the reason 
why such a characteristic comes out is not 
clear. Also, in the case of a tapered helix, 
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3:Antenna length 
/:Cylinder length 


it is difficult to control dimensions during 
manufacturing, and moreover, when tape is 
used for the helix, this method can not be 
used. For these reasons, the straight antenna 
type is much more practical; because its 
characteristics are plain and vivid, as describ- 
ed before. 

Now that the characteristics have thus been 
clarified, it should be possible to obtain com- 
plete matching by selecting antenna and cy- 
linder lengths, i. e. by bringing the admittance 
to the center of the Smith chart. With this 
object in view, closer examination was~ made 
by means of a magnified admittance chart, 
and then Fig. 7 was obtained. Now, since it 
has been made clear by the study described 
in the proceding section that in the best 
combination, the position of the short plunger 
is in the proximity of 2mm, the admittance 
in the case of p=2 is sought. In the actual 
measurement result, the admittance points 
are, as already mentioned before, in disorder 


Fig. 7—Magnified chart of coupling 
point admittance. 


due to small random reflections, when 
matching is good. While, in this figure, the 
admittance at only 3 points in the frequency 
band are plotted in order to clarify the 
tendency. In this magnified Smith chart, the 
change in position of points at its central 
part as against the change in antenna length 
g shows, in some places, a_retrogressive 
curve. The tendency in general, however, 
shows a change upwards and downwards. 
The best combination K, obtained in the 
preceding section by the method of statistical 
analysis shows, as expected, good matching. 
However, the curves for cylinder lengths of 
both 2.5mm and 1.5mm, as seen in the figure, 
do not pass through the center of the chart. 

This shows the possibility of obtaining 
better matching by changing cylinder length. 
Then, as the second step, cylinder length of 
2.0mm was selected and studies were made 
in broadband matching. 


570 


6. Broadband Matching 


As described above, by measuring with a 
direct-viewing device the admittance of many 
samples, general effect of each dimension has 
been made clear. And, possibility is found 
to obtain the best one with the cylinder 
length of 2.0mm. The next consideration 
should be to measure the characteristics care- 
fully to obtain a broader frequency band 
width and to determine the final design. 
For this purpose, the following two samples 
were prepared. 


ET EEE EEESESSS 
| 


Sampie l | g 
M | 2.0mm 325 mm 
N |. 2.0 2.5 


With regard to these samples, close measure- 
ment was made of the admittance for a 
wide frequency range by making use of a 
waveguide type standing wave detector. No 
mention is made here of the results one by 
one, but it was confirmed that Sample N was 
the best, and the admittance points, as shown 
in Fig. 8, pass through the center of the 
chart where the position of the short plunger 
was adjusted. In fact, the admittance chart 
of N is remarkably good as shown in Fig. 9. 

The movable short plunger is liable to 
cause many troubles such as parasitic re- 


Fig. 8—Distribution of addmittance of sample 
M and N. 
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Fig. 9—Broadband admitance of sample N1.s. 


sonance or poor contact. Therefore, it would 
be desirable to use fixed short plungers for 
the finished tubes. The effect due to the 
position of short plunger was investigated. 
Fig. 10 shows the voltage standing wave 
ratios measured, with respect to Sample N, 
at intervals of 100 KC within the range from 
10.2Ge to 12.5Gc. If the position of the 
short plunger is selected properly in such a 
way, very good matching of the standing 
wave ratio below 1.15 may be obtained for 
a band as broad as 2000 Mc. If a standing 
wave ratio of up to 1.3 is allowed, the band 
width will cover 2200 Mc, i.e. 20% of the 
central frequency. It is also confirmed from 
this figure that, within the range of 10.7-11.7 
Ge, used for 11 Ge radio relay, the position 
of the short plunger may be at any point 
between 1.4mm and 1.85mm. _ Therefore, 
even if the position of the plunger is fixed, 
satisfactory matching may be obtained. 


7. Conclusion 


A systematic examination has been made 
with respect to the coupling circuit between 
the helix and the waveguide having the 
structure shown in Fig.1. As a result, the 
effect on the matching characteristics due to 
the change of antenna length, cylinder length 
and position of short plunger, has been clari- 
fied distinctly. 


First, by the method of statistical analysis 


— 
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Fig. 10—Measured standing wave ratios of sample N at various short plunger positions. 


of the orthogonal array layout of experiments, 
it could be deduced that no _ interaction 
occurred between the effects exerted by these 
three factors on the standing wave ratio over 
the entire bandwidth and that the matching 
of the standing wave ratio of below 1.3 
without adjustment might be - realized . with 
statistical safety. It was also shown how the 
change in dimensions of each part affects the 
averaged standing wave ratio and what com- 
bination of such dimensions give the best 
matching. It may further be concluded that 
the accuracy of each dimension of + 0.lmm 
or so will not be seen so significants and 
that it is promissing for industrial production. 

Also it was ascertained by examining on 
Smith chart the admittance at the coupling 
as seen from the waveguide that antenna 
length mainly affects:the conductance factor 
and the cylinder length mainly affects the 
susceptance factor, while the position of short 
plunger affects both. Since such characteristics 
have been clarified, compensations can easily 
be taken by measuring the admittance, even 
if the matching condition is changed due to 
partial change in dimensions or raw materials 
or to thermal distortion and the like. 

With the best design thus obtained by 


applying these results and with the position 
of the short plunger finely adjusted, very 
good matching characteristics for a bandwidth 
of 20% of the center frequency could be 
obtained. Such matching characteristics are 
considered to be the best that have so far 
been reported. It was also confirmed that 
the fixed short plunger was available with 
enough safety. This, together with the sta- 
tistical deduction made before, enables us to 
be fully convinced of a bright prospect for 
mass production of this coupling. 

Such a successful result may not be obtain- 
ed either by the statistical analysis of the 
orthogonal array layout of experiments or by 
the admittance chart method alone. If at- 
tention were paid only to admittance, there 
would be innumerable combinations of dimen- 
sions which would effect matching with re- 
spect to a certain frequency. In order, how- 
ever, to determine the best combination, the 
method of statistical analysis should be 
applied. 

At the same time, by this statistical method 
alone, the effect of each dimension on ad- 
mittance chart could not be made clear. The 
best design was thus obtained by the simul- 
taneous use of the two methods. Moreover, 
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as already mentioned before, the cause of 
the above-described admittance characteristics 
can be reasoned quite clearly, so that these 
results may be applied to designs in future, 
since they are considered to exist also in the 
cases of other frequency or different di- 
mensions of helix etc. 
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